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ABSTRACT 


An  analysis  of  the  coverage  pattern  of  transmitting  equipment 
in  the  interference  region  requires  a  solution  to  the  geometrical  ray 
configuration  which  is  valid  within  the  effective  line-of-sight  of  the 
transmitter.  The  limiting  parameter  on  the  accuracy  required  in  the 
analysis  is  quite  high  for  relatively  short  wave  lengths  of  transmitted 
signal.  The  effects  of  the  atmosphere  and  physical  irregularities  will 
be  shown  to  modify  an  assumed  geometric  configuration.  In  addition, 
the  results  of  the  analysis  will  be  shown  solvable  on  a  digital  computer. 
The  techniques  of  the  solution  will  be  indicated. 
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CHAPTER  I 


INTRODUCTION 


An  analysis  covering  the  area  of  pattern  propagation  will  deal  with  the  area  of 
coverage  of  an  airborne  platform  housing  electronic  equipment  which  transmits  and 
receives  electromagnetic  waves  of  various  types  of  modulation.  The  object  of  this  paper 
will  not  be  to  analyze  the  airborne  equipment  itself,  but  rather  to  deal  with  the  pattern 
of,  or  area  of  coverage  of,  the  equipment  in  the  earth's  atmosphere,  given  the  trans¬ 
mitting  and  receiving  field  strength  functions  of  the  platform.  In  other  words,  the 
properties  of  the  airborne  antenna  will  be  an  independent  variable  in  determining  the 
electromagnetic  field  strength  at  a  given  point  in  space  relative  to  the  airborne  system. 

In  general,  the  area  to  be  covered  will  be  an  analysis  of  electromagnetic  wave  propagation 
from  the  time  it  leaves  the  airborne  platform,  until  it  reaches  its  target,  and  in  the  case 
of  a  radar  system,  returns  to  the  platform  again.  Or  for  example,  in  the  case  of  a 

radar  platform  located  at  x  .  y  above  a  reference  to  the  earth,  at  what  coordinates  will 

o  o 

a  target  be  in  the  detection  range  of  the  radar  ? 

This  paper  will  cons'der  two  methods  by  which  energy  can  reach  a  target  from  a 
transmitting  platform  and  in  general  return  again.  One  method  is  by  direct  line-of-sight 
from  the  platform  to  the  target.  The  second  is  by  considering  a  path  which  reflects  from 
the  surface  of  the  earth.  With  both  paths  of  transmission  present  at  the  same  time, 
there  will  certainly  be  phase  cancellation  and  reinforcement  of  the  transmitted  energy  at 
the  target.  This  fact  alone  means  that  if  this  effect  is  going  to  be  taken  into  account, 
then  a  precise  geometrical  analysis  is  required  if  phase  considerations  are  to  be  pre¬ 
served  and  meaningful.  However,  the  validity  of  the  geometrical  representation  must 
be  analyzed  in  order  to  adhere  to  the  physical.  For  when  the  energy  is  of  a  wave  length 
in  the  vicinity  of  one  meter,  then  if  the  geometrical  situation  described  above  is  not 
presented  accurately  by  precise  equations,  the  end  result  becomes  entirely  invalid.  To 
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this  end,  much  time  and  effort  will  be  extended.  Typical  properties  which  must  be  taken 
into  account  are: 

( a)  Divergence  of  energy  from  a  spherical  reflector. 

(  b)  Effect  of  the  radius  of  the  earth. 

(c)  Atmospheric  variations  in  the  area  involved. 

( d)  The  effect  of  polarization  on  reflection. 

Other  properties  which  will  be  analyzed  are: 

(a)  Reflector  irregularities. 

(  b)  Ray  tracing  validity. 

The  author  will  perform  an  analysis  of  the  above  in  relation  to  field  strength  versus 
spacial  position  from  the  platform  by  compiling  equations  which  represent  the  factors 
involved.  While  this  may  appear  as  a  mere  task  on  the  surface,  if  one  looks  further  into 
the  problem,  he  soon  finds  himself  with  a  set  of  equations  which  become  extremely  cum¬ 
bersome  to  solve  without  making  simplifications  which  in  turn  render  the  results  invalid 
insofar  as  phase  considerations  are  concerned.  Hence,  the  author's  approach  in  de¬ 
riving  valid  equations  which  are  solvable,  within  the  accuracy  required,  on  a  digital 
computer  such  as  the  IBM  7090  will  be  covered  in  the  paper.  The  equations  finally  used 
will  be  discussed,  as  well  as  other  possible  approaches  which  were  rejected  for  one 
reason  or  another.  Of  considerable  practical  interest  will  be  the  resulting  general  com¬ 
puter  program  to  solve  the  composite  equations. 
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CHAPTER  II 

FUNDAMENTAL  CONCEPTS  OF  PROPAGATION 

An  analysis  of  any  airborne  equipment  will  usually  involve  the  properties  which 
will  be  discussed  in  the  pages  that  follow.  Consider  first  the  transmission  of  electro¬ 
magnetic  energy  in  a  theoretical  free  space.  This  would  consist  of  transmission  in  a 
vacuum  medium  between  objects  j  remote  from  other  influences  that  transmission  is 
essentially  unaffected. 

A.  Generalized  Philosophy 

In  general,  first  to  be  considered  is  the  isolated  radiating  antenna  itself.  Consider 
a  polar  axis  of  a  spherical  coordinate  system  such  that  its  origin  is  the  origin  of  the 
antenna.  A  point  in  space  relative  to  the  antenna  is  then  characterized  by  specifying 
angles  0  and  <p  with  respect  to  the  antenna.  9  may  be  the  zenith  angle  and  <p  the  angle 
about  the  polar  axis. 

An  antenna  has  radiating  properties  which  are  thus  specified  in  terms  of  9  and  <p  . 

F (9,  <p)  may  then  be  defined  as  the  ratio  of  electric  or  magnetic  field  strength  in  the 
(0,  <p)  direction  to  that  which  occurs  in  the  direction  of  maximum  field  strength.  Meas¬ 
urements  are  taken  at  a  distance  which  is  large  compared  to  both  the  wave  length  and 
dimension  of  the  antenna  such  that  the  near  field  of  the  antenna  is  neglected  and  the 
energy  transmitted  consists  of  constant  phase  wave  fronts. 

Considering  radiation  from  the  antenna,  Poynting's  vector  (P  =  E  x  H).  or  the 
cross  product  of  the  electric  and  magnetic  field,  gives  the  magnitude  of  the  energy  flow 
per  unit  area  in  the  direction  of  propagation  when  its  magnitude  is  averaged  over  time. 
For  the  steady  state  sinusoidal  consideration  this  average  power  can  be  expressed  in 
complex  notation  as : 

P  =  £  Re(E  x  H*) 
av  2  v  ' 


(1) 
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The  above  notation  assumes  that  the  incremental  section  of  the  radiation  can  be  con¬ 
sidered  as  a  plane  wave  front,  which  is  true  in  the  far  field  of  the  antenna.  Utilizing 
MKS  units,  the  instantaneous  value  of  Poynting's  vector  is  given  as: 


P  = 


(2) 


where  E  and  H  are  the  instantaneous  values  of  the  total  electric  and  magnetic  field 
strength  and: 


•V5 


ai  120ff  Ohms 

impedance  of  free  space 

(3) 

-7 

=  4!T  X  10 

henrys  /  meter 

(4) 

_  10 '9 
e0  36  Jr 

farads/ meter 

(5) 

Mq  is  called  the  permeability  and  t ^  the  permittivity.  The  time  average  of  Poynting's 
vector  hence  degenerates  to : 


as: 


P 


av 


(6) 


F(0 ,  0) 


E(0.<fr) 

E„ 


(7) 


V*'*1'  lF("',>|25av(0> 


(8) 


where  Pav  is  the  average  power  in  the  direction  of  maximum  transmission. 


At  this  point,  it  is  useful  to  define  an  isotropic  radiator  as  an  antenna  which  radiates 
equally  in  all  directions,  hence  for  all  values  of  8  and  <p: 

F(0,  0)  =  1. 0  (9) 

For  an  isotropic  radiator,  Poynting 's  vector  is  given  by  the  well-known  relationship: 


P 

iso 


W. 

_ i_ 

2 

4  ttR 


(10) 
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where  Is  the  average  power  radiated  by  the  isotropic  element  in  watts  and  R  is  the 
distance  in  question. 

The  power  radiated  by  a  directional  antenna  can  be  obtained  from  Equation  (8). 
Denote  this  power  as  W  ’  .  The  power  in  question  is  found  as  the  surface  integral  of 

A 

Equation  (8)  or  more  specifically  as: 

w2=  j*  Pav(8.  0)R2dcr  (11) 

4ir 

where  the  surface  in  question  is  a  sphere  of  radius  R  from  the  radiating  element  with 
2 

surface  area  R  dc  where  <r  is  the  solid  angle  in  steradians  as  viewed  f-^m  the 
antenna. 

With  the  above  in  mind,  the  absolute  antenna  gain  can  be  defined  and  obtained.  The 
absolute  antenna  gain  is  defined  as  the  ratio  of  the  radiated  power  by  the  isotropic 
element  required  to  produce  a  given  value  of  Poynting's  vector  at  a  set  distance  to  the 
power  radiated  by  the  directional  element  producing  the  same  value  of  Poynting's  vector 
at  the  same  distance  in  the  direction  of  maximum  power  radiation.  Rewriting  Equation 
(11)  for  a  directional  element  yields : 


W2  ~  4'R\v(0)  I4Jr 
The  absolute  antenna  gain  is  hence  equal  to: 

P  \  d<r 
W.  iso  J 

_ i_  _  _ 4v _ 

Gt  "  W2  "  r  2 

da 


(12) 


(13) 


By  employing  the  above  definition,  P^maxj  of  the  directional  antenna  is  equal  to  PigQ 
and  hence  W  =  GW.  This  leads  to  the  conclusions  that  the  maximum  values  of 

1  t  Lt 

Poynting's  vector  for  the  directional  case  is  equal  to: 


P 

max 


W. 


4m 


W2°t 

47TR2 


(14) 
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As  W  Is  the  transmitted  power,  the  substitution  that  W  =  W  can  be  made.  This  leads 
2  2  t 

to  an  expression  for  P(0 ,  0)  by  utilizing  Equation  (8)  and  making  the  substitution  that 

P  =  P 

(max)  av(0) 

W  G  2 

P  ( 6 ,  0)  =  I  F(0,  0)  |  (15) 

av  4*R2 

Through  the  use  of  the  cross  product  representation  for  Poynting 's  vector,  Equation  (2), 
values  of  the  electric  field  strength  as  a  function  of  position  are  found  as : 

-  /60W^ 

E(0,  0)  =  V  Pav120ff  = - - -  F(0,  0) 


E(fl,  0)  =  E0F(e,  0)  (16) 

It  should  be  noted  that  if  the  Hertzian  dipole  or  the  half-wave  doublet  are  to  be  considered 
as  the  gain  standards,  then  the  values  of  total  gain  should  be  divided  by  3/2  and  1.64 
respectively. 


B-  One-Way  Transmission 

For  the  preceding  information  to  be  of  practical  use,  the  transmitted  energy  must 
be  received  by  a  receiving  element  and  receiver.  If  we  consider  a  properly  matched 
receiver  and  receiving  antenna  which  are  oriented  in  a  direction  for  maximum  reception, 
the  received  power  W  can  be  written  as  ^ : 


W  =  P  A 


r  av  e 

where  P  is  the  value  of  the  received  Poynting 's  vector  and  A  the  effective  cross 
av  e 

section  of  the  receiving  antenna.  With  the  wave  length  of  the  transmitted  energy 
expressed  in  meters  (A),  the  gain  of  the  receiving  element  is  given  as: 


(17) 


4?rA 


(18) 


[1]  S.  A.  Schelkunoff,  Electromagnetic  Waves,  Van  No  strand,  New  York,  1943,  Sec.  9.4. 
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If  the  receiving  and  transmitting  antennas  are  of  the  same  polarization  and  if  the 
direction  of  the  receiving  array  is  specified  relative  to  the  incoming  energy  in  terms  of 
$'  and  <p\  then  the  received  power  obtained  by  Equation  (17)  as  a  function  of  6'  and  <p' 
must  be  multiplied  by  the  square  of  the  antenna  pattern  function,  F(6'  ,<P’).  This  leads  to 
the  following : 

P  G  X2  2 

wf=  av-7r— -  If (*•,*•)!  (i9) 

or: 

\V  G  PG  A2  2  2 

^  „  IF  («'.*')  I  I  F(0,  «#>)  | 

t  P4  ir4ifBT 


W 
_ r 

W, 


G.G  A 
t  r 


F(fl,  <(>)  F(0 ',  <f>') 


(20) 


t  (4irR) 

For  the  case  stated  previously  where  both  antennas  are  aligned  for  maximum  conditions, 
Equation  (20)  degenerates  to : 


W 
.  1 
W. 


t 


Solving  for  R  in  Equation  (21)  results  in: 


G.G  A 
t  r 

(4?rR)2 


(21) 


R 


I  w.  /gg 

_  /  _ t_  V  t  r 

W  4 it 


(22) 


This  equation  may  be  converted  directly  to  the  maximum  free-space  one-way  trans¬ 
mission  range  or  the  target  range  at  which  a  usable  energy  is  received  by  the  receiver 
by  letting  G^  =  and  utilizing  the  minimum  useful  receiver  power.  Equation  (22) 

hence  is  transformed  into: 


K  .[51 G5 

0  ’  Wmln  4’ 


t  max 


(23) 


where  RQ  is  the  maximum  theoretical  range.  This  equation  is  used  frequently  as  a 
standard  of  performance. 
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C.  Radar  Transmission 


As  this  paper  will  be  ultimately  concerned  with  two-way  transmission  as  well  as 
one-way  transmission,  a  match  for  Equation  (23)  will  now  be  derived.  Two-way  trans¬ 
mission  infers  the  concept  of  reflecting  transmitted  energy  by  a  remote  object  from  the 
transmitter  and  receiving  a  portion  of  this  reflected  energy  at  the  transmitter.  This 
type  of  transmissionor  radar  transmission  utilizes  the  properties  of  the  radar  target. 

It  is  first  necessary  in  the  development  of  radar  transmission  to  specify  the  radar  cross 
section  or  back-scattering  cross  sectional  area.  Let  cr  be  the  dimension  of  the  back- 
scattering  area.  •  It  will  not  be  the  object  of  this  paper  to  specify  either  a  mathematical 
or  experimental  determination  of  cr  but  only  to  assume  that  its  fictional  mathematical 

value  is  known.  A  definition  of  a  is  the  area  intercepting  that  amount  of  power  which, 

r  2] 

when  isotropically  scattered,  produces  an  echo  equal  to  that  observed  from  the  target  1 


Utilizing  the  previous  definition,  the  value  of  the  back-scattered  Poynting's  vector 
at  a  point  from  the  target  is  related  to  a  as  well  as  the  incident  Poynting's  vector  on  the 
target  in  the  following  manner: 


'  bs 


V 

4TTR2 


or: 


(24) 


cr 


P. 


(25) 


l 

where  P.  is  the  incident  value  and  P^g  the  back- scattered  value  of  the  Poynting's  vector 
respectively.  R  is  the  distance  from  the  radar  target.  In  a  manner  similar  to  that 
previously  discussed,  the  received  power  by  the  radar  antenna  is: 


W 


Fr(0',<*>) 


A  P.  cr 


4rrR 


(26) 


[2]  D.  E.  Kerr,  Propagation  of  Short  Radio  Waves,  McGraw-Hill,  New  York. 
1951,  p.  33. 
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The  value  of  the  incident  Poynting's  vector  is: 

W.G  2 

P,  =  — ,  I  F(0.  <P)  I 

1  4irR 2 

Assuming  that  the  transmitting  and  receiving  antennas  for  radar  transmission  are 
identical: 


W 

=  lFr(«.0) 


2  2 
I  |F(fl.0)| 


A  cr  G 
e  t 

2  2 
4*R  47rR 


where: 


Substituting  the  value  for  Ag  yields: 


W 


r 


G2X2a 

(4<r)V 


4 

I  F(®.  0)  | 


(27) 


(28) 


As  before,  afree-space  radar  range  is  defined  as  that  range  which  produces  a  minimum 
radar  signal  which  is  useful.  Solving  Equation  (28)  for  this  range: 


Wt  G2*2. 


W 


I  F(0.  <p) 


m 


(29) 


» JUL 

0  Y  Wr4  7T 

Utilizing  the  preceding,  a  free-space  coverage  factor  may  be  determined. 
Equation  (30)  into  Equation  (29)  yields  : 


(30) 

Substituting 


R  =  R0  |F(0,  0)|  (31) 

This  expression  in  general  holds  for  both  one-way  and  radar  transmission.  Utilizing 
the  coverage  factor,  afree-space  coverage  diagram  may  be  determined.  This  coverage 
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diagram  encloses  a  volume  Inside  which  the  electromagnetic  field  Is  larger  than  the 
smallest  useful  value. 

D.  The  Propagation  Factor 

Up  to  this  point,  the  free-space  condition  has  been  the  major  consideration.  In  a 
practical  system,  the  coverage  of  transmitted  energy  may  be  largely  affected  by  sur¬ 
rounding  physical  conditions.  If  reasonable  parameters  are  used  in  the  equations  above 
and  then  the  actual  ranges  determined  from  experimentation,  there  will  be  considerable 
discrepancy  in  the  magnitudes  of  the  quantities  involved.  The  main  reason  for  this 
discrepancy  is  that  operation  in  the  atmosphere  close  to  the  earth  does  not  compare  with 
the  conditions  assumed  for  free-space  transmission.  In  an  effort  to  analyze  the  true 
situation,  the  free-space  relationships  derived  will  not  be  disregarded  but  instead  they 
will  be  used  as  the  building  block  in  describing  the  new  or  actual  situation. 

Consider  first  a  factor  which  shall  be  called  the  pattern-propagation  factor.  This 
factor  is  defined  as  the  ratio  of  the  amplitude  of  the  electric  field  at  a  given  point  under 
certain  conditions  to  the  amplitude  of  the  electric  field  under  the  free-space  ideal  condi¬ 
tions  where  the  beam  of  the  transmitter  is  directed  to  the  point  under  consideration.  In 
equation  form,  the  above  is  restated  as: 


F  = 


E 

E„ 


(32) 


where  F  is  the  propagation  factor,  E  the  electric  field  at  the  point  in  question,  and  EQ 
the  reference  free-space  electric  field  having  the  same  polarization  as  E.  F  can  be  seen 
to  be  entirely  independent  of  internal  properties  of  the  transmitter  and  dependent  only 
upon  the  properties  of  the  path  of  transmission  and  the  directive  properties  of  the  antenna. 

In  the  presence  of  an  earth,  F  is  the  effect  of  the  following:  (1)  the  diffraction 
phenomenon  or  where  the  earth's  shadow  gives  rise  to  diffraction  of  electromagnetic 
energy;  (2)  the  refraction  phenomenon  or  the  effects  of  an  inhomogeneous  atmosphere 
about  the  earth,  and  (3)  the  interference  phenomenon  which  is  the  effect  of  the  earth  re¬ 
flecting  and  scattering  radiation  which  produces  an  interference  pattern.  The  fourth 
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factor  which  F  takes  into  account  is  the  directional  characteristics  of  the  antenna;  or 
for  a  free-space  consideration  F  would  be: 

F  =  !  F(0,  <*>)  |  (33) 

Hence,  the  one-way  transmission  equation  may  now  be  written  in  terms  of  the  above 
effects  as: 


W  GG  X  2 

t  (4ffR) 


(34) 


Notable  is  the  fact  that  F2  is  the  ratio  of  the  magnitudes  of  the  Poynting's  vector  at  the 
point  in  question  to  the  magnitude  of  the  Poynting's  vector  whose  maximum  free-space 
value  corresponds  to  the  same  absolute  distance  as  the  point  in  question. 


Not  quite  so  easy  to  verify  is  the  application  of  F  into  the  radar  range  equation  in 

order  to  make  calculations  that  are  valid  for  regions  other  than  free  space.  Using  the 

1 3] 

theory  applied  by  H.  A.  Lorentz  in  the  development  of  the  reciprocity  theorem  it  can 
be  shown  that  the  isotropic  source  equivalent  to  the  radar  target  is  (W^G^/  4ttR  )  F  . 
Also,  using  this  theorem  it  can  be  proven  that  the  propagation  factor  squared  is  the 
quantity  that  converts  the  free-space  transmitted  energy  into  the  energy  that  would  come 
from  the  target,  if  the  target -reflected  energy  were  incident  on  the  maximum  antenna 
pattern  of  the  receiver.  Hence,  the  radar  transmission  equation  can  be  rewritten  as: 


W 


r 


G2A2a  4 
(47i)  3R4 


(35) 


In  Equations  (34)  and  (35)  are  the  theoretical  capabilities  for  determining  field 
strength  calculations  in  the  presence  of  the  earth  and  its  atmosphere.  These  equations 
are  subject  only  to  the  conditions  that  the  transmitted  energy  be  incident  upon  either 
antenna  or  target  as  a  plane  electromagnetic  wave. 


While  Equations  (34)  and  \,35)  look  relatively  easy  to  solve  in  nature ,  it  must  not 
be  forgotten  that  the  secret  to  their  practicality  is  a  correct  formulation  of  F.  In  general, 
it  is  to  this  desired  end  that  a  large  portion  of  this  paper  will  be  concerned. 


[3]  Ibid.  ,  pp.  693-699. 
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First  consider  Figure  1.  The  presence  of  the  earth  can  be  immediately  seen  to 
produce  scattering  of  the  transmitted  energy  from  the  transmitter  located  at  a  height  of 
Zl  above  the  earth.  This  gives  rise  to  an  interference  phenomenon  at  the  receiver  or 
target  located  at  a  height  Z  above  the  earth.  From  the  knowledge  of  electromagnetic 


Figure  1 .  Roy  Representation  of  Propagation  in  the  Interference  Region  Above  A 
Curved  Earth 
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propagation  it  is  immediately  seen  that  the  energy  arriving  at  point  (2 )  in  space  is  made 

up  of  a  direct  energy  wave  front  and  a  reflected  energy  wave  front  which,  due  to  their 

difference  in  distances  traveled,  are  not  of  the  same  phase.  Or,  in  terms  of  the  E  and 

H  vector  representation,  the  vectors  representing  the  direct  ray  have  not  made  as  many 

complete  revolutions  as  have  the  vectors  representing  the  scattered  ray.  For  a  more 

technical  discussion  of  the  theory  represented  above,  the  reader  is  referred  to  any 

[4) 

standard  text  on  field  and  wave  theory. 

Hence,  in  view  of  the  preceding,  it  can  be  seen  that  the  path  length  difference 

between  the  direct  transmitted  energy  and  the  scattered  transmitted  energy  received  by 

the  receiver  in  the  case  of  one-way  transmission, or  the  equal  resulting  effect  in  the  case 

of  a  radar  type  transmission, must  be  considered  as  a  component  part  of  F.  Or  more 

explicitly,  in  order  to  show  the  qualitative  properties  of  the  pattern-propagation  factor 

F  in  the  interference  region, consider  again  the  situation  shown  as  Figure  1.  Here  will 

be  seen  the  transmitting  antenna  having  a  vertical  plane  pattern  of  F(  0)  located  at  a 

height  Z  above  a  spherical  earth.  The  angles  0,  and  <p  are  therefore  needed  to  de- 

termine  the  antenna  patterns  f(«p  )  and  f(0„)  in  the  direction  of  the  direct  and  scattered 

1  6 

rays  respectively.  Note  that  it  has  been  assumed  that  the  energy  wave  front  can  be 
represented  as  a  pencil-like  ray  or  composed  of  vectors  normal  to  the  wave  fronts  and 
in  the  direction  of  propagation.  The  validity  of  the  assumption  will  be  discussed  at  a 
later  time. 

Consider  that  in  the  direction  of  maximum  propagation  of  the  antenna  the  electro¬ 
magnetic  field  due  to  the  direct  ray  is  E^.  Then  it  follows  that  from  the  definition  of 
f(0)  =  F( 0.  0)  that  the  value  of  the  direct  ray  at  an  angle  O ^  (where  the  horizontal 
dependence  has  been  omitted)  is  : 

Ed  =  V«V  (36) 

where  and  <p^  are  relative  to  the  normal  of  Z^.  The  distance  traveled  try  the  direct 
ray  to  the  receiving  antenna  is  H0 . 


[41  S.  Ramo  andJ.R.  Whinnery,  Fields  and  Waves  in  Modern  Radio,  John  Wiley  and 
Sons,  New  York,  1953,  pp.  270-314. 
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In  a  similar  manner,  the  value  of  the  reflected  wave  as  it  leaves  the  transmitting 
antenna  is : 


Er  V(02)  (37) 

The  path  length  of  the  reflected  or  scattered  ray  is  R  +  R  .  The  direct  and  scattered 

i  O 

rays  thus  travel  a  path  length  difference  of: 

AR  =  Ri  +  R3-R2  (38) 

where  AR  >  0  for  all  locations  of  the  receiver  in  the  interference  region.  Because 
of  this  path  length  difference,  the  scattered  ray  is  retarded  in  phase  by  a  value  (AR)A 
where  A  is  the  quantity  which  converts  the  path  length  difference  into  its  equivalent 
difference  in  terms  of  the  wave  length  of  the  propagated  energy. 

The  reflected  wave  also  undergoes  a  phase  retardation  because  it  is  reflecting  from 
the  surface  of  the  earth.  Not  only  does  the  reflected  wave  suffer  retardation  but  its 
magnitude  is  affected  by  the  reflection.  At  this  point,  define  T  as  the  ratio  of  the  inci¬ 
dent  field,  electric  or  magnetic,  to  the  reflected  field.  This  is  a  complex  quantity 
expressing  both  the  phase  and  magnitude  changes  suffered  by  a  wave  of  energy  reflecting 
from  the  surface  of  the  earth.  Or,  in  terms  of  symbolic  complex  notation: 


where  p  expresses  the  fraction  of  the  incident  field  that  is  reflected  and  <p  is  the  angle 
by  which  the  incident  energy  is  retarded. 

If  the  difference  in  path  lengths  is  expressed  as  the  equivalent  phase  retardation, 
A(AR),and  4>  expressed  in  a  similar  manner,  then  the  total  phase  retardation  of  the 
reflected  complex  field  relative  to  the  direct  complex  field  is  given  as : 

A^  =  A(AR)  +  0  (39) 

where  A  tp  expresses  the  total  phase  retardation. 

If  we  take  into  account  the  divergence  of  energy  due  to  its  incidence  on  a  spherical 
reflector  and  this  quantity  is  denoted  as  D,  then  the  total  electric  field  strength  at  the 
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receiving  point  located  at  Z  above  the  earth  can  be  expressed  as : 

it 

E  =  E  +  E 
t  r  d 

Or  in  terms  of  our  previous  notation: 


(40) 


+ 


T  De 


-ia(ar) 


m2) 


(41) 


Hence,  within  its  limitations  and  assumptions,  the  ratio  of  E  to  Eq  gives  the  interference 
pattern  as  formed  in  space,  or  the  pattern  propagation  factor. 


By  limitations  and  assumptions  is  meant  a  wide  scattering  of  theoretical  and 
physical  effects  which  for  propagation  in  an  atmosphere  have  been  neglected  in  the  pre¬ 
vious  analysis.  Most  of  these  effects  will  be  discussed  in  Chapter  V  of  this  paper. 
However,  there  are  two  factors  which  will  be  mentioned  at  this  point  as  no  further 
consideration  will  be  made  of  them.  First,  at  the  receiving  pointjthe  energy  wave  fronts 
will  not  be  arriving  at  the  same  directions  as  has  been  assumed  in  the  preceding  analysis. 
The  direct  and  scattered  wave  fronts  will  in  actuality  be  traveling  in  slightly  different 
directions.  However,  as  found  from  experimentation,  this  is  important  only  when  the 
energy  wave  front  is  vertically  polarized  and  the  reflection  and  incident  angles  with 
respect  to  the  earth  differ  to  a  large  extent.  In  the  paper  these  angles  will  be  made  or 
assumed  equal  according  to  the  well  known  law: 


sin«  shut 

i-  it 


where  v  and  v  are  the  velocities  of  propagation  in  the  incident  and  reflected  mediums 

1  it 

and  a  and  a  the  incident  and  reflected  angles.  As  the  incident  and  reflected  waves 

1  it 

are  traveling  in  the  same  medium,  v  =  v  and  hence: 


«1=a2  (42) 

The  second  factor  which  will  be  neglected  is  the  small  spatial  attenuation  difference 
encountered  by  the  direct  and  reflected  energy.  In  the  large  majority  of  practical  cases, 

i  5] 

this  effect  has  been  found  to  be  negligible 


[5]  L.N.  Ridenour,  Radar  System  Engineering,  McGraw-Hill,  New  York,  1947,  pp.  58-62. 
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CHAPTER  IH 

THE  FORMULATION  OF  THE  DIFFERENCE  BETWEEN 
THE  PHASE  OF  DIRECT  AND  SCATTERED  ENERGY  DUE 
TO  PATH  LENGTH  DIFFERENCE  —  AR 


Referring  to  Equation  (41),  it  is  seen  that  one  factor  which  must  be  found  is  AR. 
Musi  uf  the  remainder  of  this  paper  will  be  concerned  with  finding  these  variable  para¬ 
meters  as  a  function  of  ZQ,  Z^  and  R^ .  The  reason  for  this  is  that  the  desired  end 

result  is  a  coverage  diagram  represented  in  terms  of  Z  and  R  ,  with  Z  as  a  fixed 

&  Z  1 

parameter.  There  has  been  considerable  effort  devoted  by  various  concerns  and 
individuals  in  the  determining  of  A  R  as  a  function  of  the  previously  mentioned  variables. 
Some  of  these  methods  will  now  be  presented  as  a  background  for  the  final  method  used 
in  this  paper. 


A.  The  Exact  Equations 

The  geometrical  situation  described  in  Figure  1  on  page  12  can  be  represented  in 
exact  forms  by  the  following  set  of  equations: 

a1=«2  (43) 

/z.  \ 

E^-l+|  —  )i  -  1,  2  (an  intervening  parameter)  (44) 


\  e/ 

R2  =  Ae  J  E  2  +  E22  -  2EjE2  cos  6l  (45) 

4E  E  cosV  -(4E  E2  1  4E2E  cos  0  )cos\  +  (E2  +  E2  +  2E  E  cos  0  -  4E2E2  )cos2  0 

1  u  m  1  w  1  Z  I  Z  I  A  1  Z  J.  1  m  Z 

M4E  E^cosV  f  2E  E,2sin2  0  *■  4E2E  cos  0,)cos  0  (46) 

i  u  1  1  u  1  1  A  1  Z 

1  (E^E2sin“  01  -  E2  cos2  9  -  E2  -  2E:E2cos  OJ  =  0 

Z ’  =  (A  f  Z  )cos  0,  -  A 
1  e  1  2  e 


(47) 
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(48) 

(49) 

(50) 

(51) 

(52) 

(53) 

(54) 

(55) 

(56) 

(57) 

(58) 
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These  equations  are  referred  to  as  exact  in  that  they  represent  the  geometrical 
situation  represented  by  Figure  1.  Theoretically  they  may  not  be  exact  because  they 
deviate  from  what  is  physically  happening  or  because  they  cannot  be  solved  to  the  re¬ 
quired  precision.  The  derivations  of  all  formulae  are  straightforward  from  the 
geometrical  viewpoint  and  may  be  checked  fcy  the  reader  if  desired. 

The  use  and  practical  application  of  these  equations  may  be  questioned.  For  an 
explicit  solution  of  A  R  as  a  function  of  Z  ,  Z  and  R  one  would  first  proceed  to  solve 

X  Z  Z 

the  quadratic  in  cos  9  from  Equation  (46) .  This  would  give  the  cos  0  as  a  function  of 
Z  2 

Z  ,  Z  ,  and  cos  9  .  Equation  (45)  would  then  be  solved  for  cos  9  .  Making  appro- 
12  1  1 

priate  substitutions,  AR  as  a  function  Z^,  Z^  and  R2  can  then  be  determined.  Equation 

(57)  gives  a  measure  of  the  section  in  which  the  equations  are  valid,  or  Z  is  in  the 

z 

interference  region.  This  means  that  if  R  is  greater  than  L,  then  Z  has  a  minimum 

z  z 

value  greater  than  zero.  While  theoretically  the  preceding  is  possible-,  at  its  best  it  is 
very  complicated.  Also,  in  order  to  maintain  the  values  R  and  R  to  the  necessary 

X  o 

degree  of  precision  in  terms  of  the  wave  length  of  propagation  involved,  the  number  of 
accurate  digits  one  is  required  to  maintain  can  be  seen  to  be  unwieldy. 

Variations  in  the  exact  geometrical  representation  or  equations  have  been  derived 
by  various  writers.  Most  of  these  involve  implicit  relations  which  indeed  are  easier  to 
work  with  than  those  described  previously  but  their  use  is  limited  in  the  plotting  of  a 
coverage  diagram  in  terms  of  the  parameters  previously  stated.  The  University  of 
Michigan^  and  the  Fairchild  Engine  and  Airplane  Corporation^  have  both  contributed 
concepts  not  generally  used  in  the  preceding  equations.  Equation  (12)  in  the  Fairchild 
report  is  questionable,  should  reference  be  made  to  this  paper.  The  author  has  not, 
however,  proved  it  incorrect  even  though  dimensional  analysis  indicates  an  error. 


[6]  T.B.  Curtz,  M.  L.  Barasch,  et  al,  "Analysis  of  Padar  and  Its  Modifications  — 
Final  Report,  "No.  2476-1-F,  April,  1956. 

17)  "Padar  Investigation  —  Final  Report, "  No.  64R-10,  June,  1956. 
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B.  The  Cubic  in 

Admittedly,  there  are  simplifications  that  can  be  made  in  the  so-called  "exact" 
equations.  One  of  these  may  be  termed  as  a  cubic  in  r^ .  Consider  the  following 
approximation: 

Z  and  Z  «  A  (59) 

1  1*  6 


Therefore: 


r,  and  r„  «  A 
1  2  e 


*  Z2 


(60) 


(61) 


(62) 


dj  *  r1’  d2  “  r2’  d“  r 


(63) 


ip  s  tan 


<2,  *  Z2> 


2  ,  2 
ri  r2 

2Ae<Zl  +  Z2>J 


r  -  2A  Z, 

r  *  1 _ S_L  +  , 

2  “  2r 


r  -  2A  Z, 
1  el 

2r. 


+  2  A  Z, 
e  2 


(64) 


(65) 


2r^  -  3rrZ  +  (r  -  2A  Z  -  2A  Z„)r  +  2A  Z  r  2  0  (66) 

1  1'  e  1  e  2'  1  e  1 

Equations  (61)  and  (62)  follow  directly  from  Equations  (47)  and  (48)  by  letting 
02 

cos  6  a  1  -  .  Likewise,  Equation  (63)  results  from  Equations  (58)  and  (49)  by 

letting  sin  6  a  0.  Equation  (64)  results  from  Equation  (50)  utilizing  Equations  (61),  (62), 
and  (63).  Equations  (65)  and  (66)  follow  from  (64)  assuming  that  (64)  is  valid.  Equation 
(66)  has  the  solution: 
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r  H±UH 


rrr  PC0Bi  3 


(67) 


where: 


P  = 


2/  Ae(Zl  +  Z2)  ♦  (§) 


(68) 


<p  =  cos 


-li  2Ae(Z2  -  Zl>r  / 


I 


(69) 


Z1^Z2 


For  Z  <  Z  a  similar  equation  can  be  found  in  terms  of  r  .  This  solution  is  essentially 

1  2  r  gj  “ 

that  given  by  Kerr  1  on  page  113  except  for  the  sign  change  in: 


r  =£  .Pcosji-^-3-! 
12  I  3  ) 


(70) 


This  can  be  shown  to  be  incorrect  if  we  observe  that  as  Z  £  Z  ,  the  following 

1  2 

inequalities  must  hold: 


ri  ~  2  ~ 


However,  as: 


cos  <p  so . 


o  o 

270  >  <p>  90 


O  4>  +  7T  O 

150  >  —  >90 


(71) 

(72) 

(73) 


COS 


(0 +  T) 


s  o 


(74) 


[8]  Kerr,  op.  cit.  ,  p.  113. 
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Hence  for  (71)  to  hold  true,  Equation  (70)  must  be  incorrect  and  the  form  must  be  that 
of  Equation  (67).  The  law  of  cosines  may  be  invoked  to  find  R  and  R  . 

1  O 

R1  =  A!  +  (Ae  +  Zl)2  "  2Ae(Ae  +  Zl>  008  °1  (75) 


e 


i 


(76) 


R2  =  A2  +  (A  +  Z  J2  -  2A  (A  +  Z  )cos  6 
3  g  e  2  g  g  2 


(77) 


e  = 


(78) 


Therefore : 


ar  =  -r2  + 


A2  +  (A  +  Z,)2  -  2A  (A  +  Z  lcos 
e  e  1  e'  e  1' 


+ 


Ae*<VZ/-2W 


Z2)cos 


(79) 


It  should  be  noted  that  the  cubic  in  r  does  not  represent  a  full  third-order  approxi¬ 
mation  to  r^.  To  obtain  this,  the  third-order  term  in  the  expansion  of  sin  0  or  sin  6^ 
must  be  used  in  the  approximation  of  d^  and  d^ .  This  will  not  only  change  the  coeffi¬ 
cient  of  r  but  will  also  change  the  remaining  coefficients.  However,  as  Z^/  A^  and 

Z  / A  are  small,  these  changes  are  unimportant.  Theoretically,  these  equations 
2  6 

break  down  to  the  largest  extent  at  the  extreme  conditions  or  where  Z2  is  too  large  or 
where  R2  is  at  either  extreme. 


Considerable  work  has  been  done  on  obtaining  the  solution  to  the  cubic  in  r  . 
Transformations  have  been  found  which  normalize  the  variables  in  such  a  way  that  the 
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computing  required  for  a  solution  has  teen  reduced  considerably.  The  application  of  this 

[9] 

transformation  is  recommended  and  may  be  found  in  Fishback  and  Burrows  and 
Atwood1101  . 


C.  The  Flat  Earth  Approximation 

A  third  method  of  attack  in  determining  A  R  in  the  interference  region  can  be 
found  by  assuming  a  flat  earth  or  an  earth  whose  radius  is  infinite.  Under  this 
assumption: 


-1 


\  =  tan 


Z  + 
1 


If  Z  /r  and  Z  /r  are  assumed  small,  then: 

11  u  it 


R  »  r  + 
11 


R„ 


(82) 


(83) 


(84) 


V 

r  a  - - 

1  Zl+Z2 


(85) 


[91  W.  T.  Fishback,  "Simplified  Methods  of  Field  Intensity  Calculations  in  the 
Interference  Region,  "  Report  No.  461,  1943. 

[10]  C.R.  Burrows  and  S.  S.  Atwood,  Radio  Wave  Propagation,  Academic  Press, 
New  York,  1949. 
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Z  r 
2  1 

r  =>  - 

2-  Z, 


(86) 


AR 


2Z1Z2 


(87) 


Equation  (87)  can  be  seen  from  the  fact  that  pure  geometrical  deduction  indicates: 


r2=  V  r  +<VZ1> 


<VR3>  =  V  r  +<VZ1> 


Expanding  R  and  (R  +  R  )  by  the  binomial  theorem  and  subtracting  yields: 

u  1  O 


,  “I  *  10Z1Z2  *  322 


2r 


8r 


or: 


AR 


2Z  Z 
12 


(Assuming  Remaining  Terms  Small) 


Needless  to  say,  these  equations  break  down  everywhere  that  the  previous  equations 
did  plus  their  additional  inherent  loss  of  accuracy  brought  about  by  the  extended 
assumptions  and  simplifications.  Notable,  however,  is  that  the  cubic  in  r^  is  now  re¬ 
placed  by  (85)  and  (86).  Of  particular  interest  is  the  fact  that  these  equations  break 
down  very  rapidly  in  the  region  where  R  approaches  or  is  greater  than  L  and  Z  is 
small. 


Reference  will  be  made  in  a  later  section  of  the  paper  as  to  the  quantitative 
magnitudes  of  the  errors  involved  in  these  and  the  previous  assumptions. 
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D.  The  Author1  s  Analysis  —  The  Ex^ct  Transcendental  Equations 

Three  methods  of  determining  AR  have  so  far  been  presented.  Each  succeeding 
approach  is  simpler  in  terms  of  the  computations  required  to  obtain  a  result.  However, 
the  simplifications  are  obtained  at  an  increasing  loss  of  accuracy.  It  is  the  author's 
intent  in  this  section  to  derive  equations  which  will  yield  the  desired  AR  to  a  high 
degree  of  accuracy  comparable  to  that  obtained  theoretically  by  the  fourth  order  "exact" 
equations  but  at  the  same  time  not  to  have  the  equations  too  unwieldy. 


Figure  2.  Vector  Representation  of  Interference  Region 
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Consider  Figure  2  which  is  essentially  the  same  as  Figure  1  drawn  in  vector 
notation.  The  vectors  V  ,  V  ,  and  V  may  be  expressed  in  the  following  manner: 

1  m  O 


-  (x  -  R  cos  0)i  +  (y  -  R  sin  0)j 


(88) 


=  (R  cos  0)i  +  (R  sin  0)j 


V  =  (x  -  Rcos  0)i  +  (-  R  sin  0)j 

O  J 


(89) 

(90) 


Forming  the  following  scalar  or  dot  products  results  in: 


Vj  •  -  xRcos  0  -  R2cos2  0  +  Rysin  8  -  R2sin2 


2  2  2 

xRcos  0  +  yR  sin  0  -  R  (sin  0  +  cos  0) 


(91) 


=  xR  cos  0  +  yR  sin  0  -  R 


V  •  V  =  x  Rcos  0  -  R2  cos2  0  -  R2sin2  8 

O  u  *} 


2  2  2 

x^Rcos0-R  (sin  0  +  cos  0) 


(92) 


=  x  Rcos  0  -  R 

w 


However,  by  the  definition  of  the  scalar  or  dot  product  of  two  vectors: 

V  V  |V1(  iv2  I  COSQ^ 


V  V|V3l|V2|cOSa2 


(93) 

(94) 


|v  Icoso  =  ^Ji£os  e  4;  yR sln  &  -  r 

11  IV 


(95) 


27 


But : 

|v2I=h 

|V^  |coso^  =  xcos  9  +  ysin 9  -  R 
|V  |  cos  a  =  x  cos  9  -  R 

o  Z  o 


(96) 

(97) 

(98) 


coscv. 


_  xcos  0  +  ysin  0  -  R 
1  ”  |V  I 


x  cos  0  -  R 
3 


IV, 


or: 


cos  a 


xcos  6  +  ysin  6  -  R 


1  / 


(x  -  R  cos  9)  +  (y  -  R  sin  9) 


cost* 


XgCOs  9  -  R 


(x  -  Rcos0)2  +  (RsinO)2 

O 


Using  the  following  three  identities : 


(1) 

X  = 

(R+Z)  cos  9^ t 

(2) 

y  = 

(R  +  Z^sin  9  ^  , 

(3) 

X3  = 

r  +  z2 

It  is  easily  proven  that: 


2 

cos  or  = 


(R  +  Z^K  -  2R  +  R2/(R  +  Z^K 
1  (R  +  Z  )l/K  -  2R  +  R2/(R  +  Z^K 


(1  +  Z^RJK  +  R/(R  +  Z  )K-2 
(1  +  Z  /R)l/K  +  R/(R  +  Z^K-2 


(99) 


(100) 


(101) 


(102) 


(103) 
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where  K  =  cos  (0  -  0  ) 
V 


(R  +  Z  )cos  0  -  2R  +  R2/(R  +  Z  )cos  0 
cos  a  =  - 2 - 

(R  +  Z  )/cos  6  -  2R  +  R  /(R  +  ZJcos  0 
2  2 

(1  +  Z  /R)cos  0  +  R/(R  +  Z  )cos  0-2 
__  _ 2  2 

(1  +  Z  /R)l/cos  0  +  R/(R  +  Z  )cos  0-2  *1<H 

2  2 

The  previous  equations  for  and  are  at  their  best  not  as  simple  as  might  be  desired. 

Assuming  that  0  is  known,  which  will  be  proven  at  a  later  time,  and  knowing  from  the 
*  2  2 

previous  discussion  that  cos  a  =  cos  a^,  the  previous  two  equations  can  theoretically 

be  solved  for  0,  hence  a  and  a  .  However,  the  solution  for  0  is  an  implicit  relationship. 
1  2 

It  would,  hence,  be  desirable  to  at  least  attempt  to  simplify  Equations  (103)  and  (104). 


With  this  goal  in  mind,  consider  the  following  triple  scalar  product: 


Wk=|V2!  lV3isin“2 


Vg  X  V1  ■  k  =  |v2  I  |V1  Isino^ 


Expanding  (105): 


V3XVk  = 


x  -  Rcos  0  -Rsin  0  0 

3 


Rcos  0  Rsin  0  0 

0  0  1 


=  (x  -  Rcos  0)  (Rsin  0)  +  Rcos  0(Rsin  0) 

u 

=  x  Rsin  0 
») 


(105) 

(106) 


(107) 


where  k  represents  the  unit  vector  in  the  third  dimension. 
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V2XV1 


Rcos0  R  sin  0  0 

x  -  R  cos  0  y  -  R  sin  0  0 

0  0  1 


R  cos  0(y  -  R  sin  0)  -  (x  -  Rcos  9)  (R  sin  0) 

=  Rycos  0  -  xRsin  0 
Utilizing  the  definition  of  the  triple  scalar  product: 


V2  X  V1  '  k  =  IV1  I  lV2  I  <  <sina1)  <cos  °°) 

=  |V1 1  |V2  |  sina^ 


(108) 


(109) 


V3  x  V2  ‘  k  =  |V3  I  |V2  1  (  (-sina2)  (cos  °°) 


Note: 


But: 


Rewriting  (93)  and  (94): 


=  IVgl  |V2I  (-sina2) 


W  V2XV1 


V1XV2  =  -<V2XV1> 


lv3l  lv2l- 


COSO! 

a 


IvJ  lv2l  = 


cosa^ 


(110) 
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.  -  _  /-flina 

•  •  V4xV  •  k  =  V  •  vj  - - 

3  2  3  2  \  cos  a 


_  /  sina 

V  x  V  •  k  =  V  •  V - - 

2  1  1  2\  cos  a. 


Therefore : 


tana. 


V_2XVk 

V  .  V 
1  2 


tana„ 


'V3X  V2  '  k 

V  .  V 
3  2 


Utilizing  Equations  (91),  (92),  (107),  and  (108): 


tana  =  ^°-s  6  '  _ 

1  2 
xRcos  8  -  yRsin  8  -  R 


ycos  9  -  xsin  B 
xcos  e  +  y  sin  0  -  R 


(111) 


(112) 


(113) 


(114) 


(115) 


As : 


-  x  R  sin  8 

tana  =  - — 

2  2 

x  Rcos  8  -  R 

O 

-x  sin  8 

—  *5 

cos  8  -  R 

(R  +  Zp  cos  0  , 
(R  +  Z^)  sin  $  , 


(116) 


and  using  the  trigonometric  identities  for  the  sine  and  cosine  of  the  sum  of  two  angles, 
the  following  is  obtained: 
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sin(01  -  0) 

tan“l  cos(S1  -  6)  -  R/(R  +  Zx) 
sin  9 

tana  -  - 

2  -cos  6  +  R/(R  +  Z  ) 

6 


(117) 


(118) 


As  transcendental  equations,  these  equations  are  relatively  simple  in  view  of  the  pre¬ 
vious  results.  From  the  previous  derivations,  the  absolute  value  of  V  (or  R^)  and 
the  absolute  value  of  V  (or  R  )  can  be  determined.  From  Equations  (97)  and  (98) : 

O 


Iv 


x  cos  6  +  y  sin  9  -  R 
cos  a^ 


(119) 


x„ cos  6  -  R 

I v  |  =  - - 

3  cos  a2 

Substituting  in  the  value  for  x,  y  and  x^ : 


(120) 


IV. 


(R  +  Z^Jcos  9  cos  9  +  (R  +  Z^)  sin  9  ^  sin  6  -  R 


cosa. 


Let 


lvn 


(R  +  Z^)  cos  9  -  R 
cosa2 


Z  =  cos  (9  -  6) 

R  =  A 

e 
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Then: 


lvi'  "Ri 


[(Z  -  1)  *  (Z;  x  Z)/Ae)]  I 

cosa^  j 


(121) 


V3  R3 


[(cos  0  -  1)  +  (Z2cos  0)  /  Ag  J  } 


cos  a. 


(122) 


To  complete  the  analysis  of  AR,  it  is  now  necessary  to  derive  an  expression  for 
6  ■  This  is  a  relatively  simple  procedure  as  it  results  directly  from  the  solution  of 
an  oblique  triangle.  Given  the  following  triangle  denoted  as  Figure  3: 


Figure  3.  Triangle  Solution  of  0^  in  General  Terms 
the  angle  included  between  two  known  sides  b  +  c  may  be  found  as: 

a  =  2tan  1  (  — - — 'l 
\  s  -  a  / 

where : 

s  =  2  (a+b  +  c) 

r=  J  (s  -  a)(s  -  b)(s  -  c) 


s 
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For  the  situation  pictured  in  Figure  1  : 


s  - 


2A  -  R  +  Z, 
e  2  1 

2 


(123) 


(8 


R2)(8 


ZjHs 


A  - 
e 


Z2> 


(124) 


6  =  2tan 


-1  <8  -  Ae  -  Zi)(s  -  Ae  -  Z2) 


(G  -  R2)S 


(125) 


The  logical  procedure  to  follow  when  solving  for  A  R  from  these  equations  is  to  first 
solve  Equations  (123),  (124)  and  (125)  knowing  A  ,  R  ,  Z  and  Z  .  Having  once  solved 

6  u  1  u 

for  6  ,  then  Equations  (117)  and  (118)  may  be  solved  for  a  and  a  .  Equation  (121) 

1  Xu 

and  Equation  (122)  may  then  be  solved  directly  for  R^  and  R^.  AR  is  then  found  as: 

A  R  =  R  +  R  -  R 
13  2 

The  only  difficulty  of  any  consequence  when  solving  these  equations  when  compared 
with  the  previous  exact  formulations  lies  in  the  solution  of  Equations  (117)  and  (118).  It 
has  been  found  by  the  author  that  these  equations  are  solved  rapidly  by  using  the  following 
technique.  First,  calculate  a  value  of  $  by  using  the  flat  earth  approximation: 


tan  ^ 


(126) 


by  the  approximation  r  a  R^ 


(127) 


(128) 
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Having  thus  found  an  approximate  value  for  9 ,  calculate  the  tangents  of  a  ,  and  by 
Equations  (117)  and  (118).  Then  find  the  difference  between  the  tangents  of  and 
and  call  this  D.  Assuming  that  D  *  0,  then  the  following  may  be  stated,  assuming  a 
flat  earth : 


(129) 


Assume  that  r/  is  the  factor  which  is  correct  and  that  0(Ag)  (l/Z1  +  l/ Z^)  is 
causing  D  to  deviate  from  zero.  The  following  then  can  be  said: 


0  A 
0  e 


+ 


(130) 


where  0Q  is  the  corrected  value  of  6.  Subtracting  Equation  (130)  from  Equation  (129) 
yields: 


V 


e 


D  +  0A 

e 


1 


or: 


0 


0 


(131) 


From  Equation  (131)  a  new  value  of  the  tangents  of  and  may  be  determined.  If 
the  difference  of  the  tangents  of  these  new  angles  differs  from  zero  within  the  accuracy 
lieing  carried  out,  then  the  author  suggests  the  use  of  the  following  technique. 


First,  assume  a  straight  line  variation  for  the  difference  of  the  tangents  of 
and  at  a  proximity  of  the  reflecting  region  of  interest.  Refer  to  Figure  4. 
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ASSUMED  IN 
THIS  DIRECTION 

Figure  4.  Straight  Line  Approximation  of  the  Difference  in  the  Tangents 
of  and  dj  vs.  the  Angle  g  . 

In  essence,  the  previous  two  equations  have  given  the  two  noted  points  0  A  and  6  A  . 

Li  Lt  11 

From  theory  stated  previously,  it  is  desired  to  have  A  0  for  absolute  accuracy. 

Hence,  where  M  is  the  slope  of  the  line  and  B  the  0  intercept: 


0  =  M0  +  B 


(132) 


where 


M  = 


<V*2> 

<vv 


b=a2  -M02 

Solving  for  the  value  of  0  0 ^  which  will  make  -  0  yields: 


A2<Vei> 

A  -  A 

2  1 


+  0. 


where  in  general: 


A2  =  (tano^  -  tanu2)  for  the  flat  earth  approximation. 


(133) 
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A  =  (tana  -  tana  )  for  the  Incremented  flat  earth  approximation  or  from 
1  1  2  Equation  (131). 


0  ^  =  the  value  of  0  obtained  from  the  flat  earth  approximation. 


0 


2 


=  the  value  of  0  obtained  from  the  incremented  fiat  earth  approximation  or 
from  Equation  (131). 


The  above  procedure  may  be  repeated  until  the  desired  degree  of  accuracy  is  obtained 

by  making  the  appropriate  substitutions  for  A^  ,  A^,  0  ^  and  0 ^  for  each  pass  through 

the  approximations.  Theoretically,  any  desired  degree  of  accuracy  may  be  obtained. 

In  practice,  the  author  has  found  the  above  procedure  to  converge  very  rapidly.  In 

general,  within  a  range  of  three  hundred  meters  and  with  altitudes  above  a  few  meters, 

the  values  of  the  two  tangents  can  be  brought  to  ten  to  twelve  decimal  places  of  accuracy 

by  no  more  than  six  repeated  approximations  of  Equation  (133).  For  values  of  Z^  and 

Z  above  two  thousand  feet  the  convergence  appears  much  more  rapidly. 

£ 

The  preceding  equations,  therefore,  provide  a  method  within  the  accuracy  of  the 
geometrical  framework,  by  which  AR  may  be  determined  to  within  any  theoretical 
desired  accuracy  ranging  from  the  flat  earth  assumption  to  the  exact  equations.  Ad¬ 
mittedly,  the  labor  involved  increases  with  the  desired  accuracy.  But,  seldom  are  the 
flat  earth  approximations  within  the  realm  of  reality  for  a  physical  situation,  especially 
when  the  transmitted  energy  in  the  interference  region  is  of  a  wavelength  in  the  vicinity 
of  ten  to  one  hundred  centimeters.  This  is  especially  true  for  small  values  of  Z^ . 

Also  obtained  has  been  an  increase  in  the  accuracy  of  A  R  for  a  specified  number  of 
significant  digits  carried  in  the  calculation.  For  all  calculations  carried  out  in  meters, 
the  analysis  of  A  R  by  the  method  outlined  under  the  heading  of  the  cubic  in  r^  will  re¬ 
quire  twenty  decimal  digits  of  accuracy  to  maintain  A  R  accurate  to  the  nearest 
one-hundredth  of  a  meter  within  the  interference  region  previously  covered.  The 
author's  approach  requires  approximately  twelve  decimal  digits  of  accuracy  to  maintain 
the  same  information.  The  author's  approach  does  not  make  any  assumptions  such  as 
are  made  in  the  evaluation  of  the  cubic  in  r^ .  hence  it  is  inherently  more  accurate. 
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CHAPTER  IV 

PROPAGATION  FACTOR  COMPONENTS  OTHER  THAN  A  R 

Having  spent  so  much  effort  in  obtaining  A  R,  which  was  previously  shown  to  be 
a  component  of  F,  one  might  reasonably  wonder  if  it  might  not  have  teen  tetter  to  find 
out  what  other  components  made  up  the  factor  F  .  If  they  proved  to  demand  a  mul'.i  - 
plicity  of  the  work  previously  described,  then  the  analysis  of  AR  might  have  been  over- 
refined.  The  overall  validity  of  the  factor  F  can  be  no  tetter  than  its  weakest  constitu¬ 
ent.  Luckily,  however,  the  determining  of  A  R  is  the  major  stumbling  block  in  finding 
the  component  parts  of  F. 

A.  The  Complex  Coefficient  of  Reflection 

Equation  (41)  indicates  that  a  second  factor  which  must  be  determined  is  F  . 

This  factor  takes  into  account  the  reflections  of  waves  from  the  earth1  s  surface.  The 
reflection  factor  was  previously  defined  as: 

p  _  -i<P 

I  =  P  t 

If  the  electromagnetic  waves  incident  to  the  reflector  are  horizontally  polarized,  then 

T  is  defined  as  the  ratio  of  reflected  to  incident  electric  field  or  the  ratio  of  the  verti- 
h 

cal  components  of  the  magnetic  field.  The  ratio  of  reflected  to  incident  horizontal 

components  of  the  magnetic  field  is  -r  .  For  vertically  polarized  waves  F^  is  defined 

as  the  ratio  of  reflected  to  incident  magnetic  field.  The  ratio  of  the  vertical  components 

of  the  electric  field  or  minus  the  ratio  of  the  horizontal  components  is  likewise  equal 

to?  . 
v 

From  these  definitions,  Fresnel's  equation  for  the  smooth  plane  earth  can  be 
derived: 


r 

v 


(KI/K0)28in(2  -  V  (K^K^'-cor  ^ 

<VKo>2  sin  i/2  +  /(K/k/  -cos2i/2 


(134) 
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-10  sin  0  -  /  (K  /  K  )2  -  cos2  4, 

T  =  P.e  h=  - ^ -  J.  0-  .. (135) 

h  h  / ^  2  2 

sin  02  +  /  (K/K0)  '  cos  02 

The  author  will  not  derive  these  equations  here  as  they  may  be  found  in  any  standard 
work  ^ 1  ^  .  A  variable  in  these  equations  is  K1  where  K3  is  defined  as  the  propagation 
factor  of  the  earth: 

-  j w/iu  (136) 

where  j  =  nT^I  ,  uj  is  the  frequency  of  the  propagating  energy,  p  the  permeability  of  the 
earth  expressed  in  henries  per  meter,  the  conductivity  of  the  earth  in  ohms  per  meter 
and  e  the  complex  dielectric  constant  of  the  earth  in  farads  per  meter.  The  properties 
of  the  earth  are  expressed  as  relative  values. 

The  variable  Kq  is  the  square  root  of  the  propagation  factor  for  air  which  for  the 
air  at  the  surface  of  the  earth  will  be  assumed  approximately  equal  to  that  of  free 
space. 

Ko=  "Vo  (137) 

where  pQ  and  £Q  are  the  permeability  and  dielectric  constants  respectively  for  free 

space.  The  permeability  of  the  earth  can  essentially  be  said  to  be  equal  to  that  of  free 

-  /  2 

space.  Hence,  in  view  of  the  above  (K^/K^)  can  be  written  as  follows: 


<R1/K0> 


c/€o  ”  j  <r/"€o 


€1*)€2 


(138) 


The  quantity  e/e^  is  the  usual  dielectric  quantity  commonly  listed  in  tables  for  dielec¬ 
tric  materials  and  ?  is  called  the  complex  dielectric  constant.  In  terms  of  the  MKS 
system  of  units: 


e2  =  o/u>co 


60 \<r 


(139) 


[11]  J.A.  Stratton,  Electromagnetic  Theory,  McGraw-Hill,  New  York,  1941,  Secs. 
9 . 4  and  9.9. 
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-  .  2 

where  X  is  the  wavelength  of  propagation  expressed  in  meters.  Hence,  ( K  /  KQ)  can 
be  simplified  to: 

7  =  -  j60Xcr 


Figure  5.  Variations  of  Complex  Coefficient  of  Reflection  With  Incident 
Angle  ^  ^  • 


(140) 


40 


In  general,  the  expected  behavior  of  Equations  (134)  and  (135)  can  best  be 
described  with  few  words  and  the  use  of  Figure  5. 

B.  The  Divergence  Factor 

A  third  factor  which  must  be  determined  is  D  as  seen  from  Equation  (41) .  This 
factor  takes  into  account  that  the  electromagnetic  energy  is  weakened  upon  reflection 
from  a  spherical  earth  by  a  divergence  of  the  concentration  of  the  energy.  The 


Figure  6.  Curved  Earth  Geometry  for  Computation  of  the  Divergence  Factor 
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following  analysis  will  help  to  clarify  the  general  concept  involved-  Reference  should 

be  made  to  Figure  6.  This  analysis  was  first  carried  out  in  detail  by  van  der  Pol  and 
[  12] 

Bremmer  in  a  detailed  analysis  of  defraction  of  waves  by  a  spherical  conductor. 

The  concept  is  based  upon  the  principle  of  comparing  the  density  of  the  reflected 
rays  from  the  earth  or  sphere  with  those  that  would  have  been  present  had  reflection 
occurred  from  a  plane  reflector.  The  field  strength  difference  is  proportional  to  the 
square  root  of  the  ray  intensity  difference.  The  analysis  is  carried  out  by  comparing 
the  cross  section  of  a  cone  obtained  from  a  spherical  reflection  with  that  obtained  from 
reflection  from  a  flat  conductor. 

The  rays  leaving  the  transmitting  antenna  have  an  effective  cross-sectional  area 
of: 

R2  sin  i L  dip.  da  (141) 

all 

where  a  is  measured  in  a  plane  perpendicular  to  the  plane  of  paper  and  R  is  the 

s 

straight  line  distance  from  the  antenna  at  Z  to  the  antenna  at  Z  .  Reflection  from  a 

1  2 

plane  earth  would  result  in  an  equivalent  cross-sectional  area  of  rays  equal  to: 

(R  +  R.)2  sin  ip.  di p  da  (142) 

This  is  because  the  reflective  path  length  appears  from  an  image  directly  below  the 
transmitting  antenna  at  a  point  -Z  .  The  cross-sectional  area  of  the  rays  in  the  spherical 
case  is  s  or  s  cos  ip  .  From  the  geometry,  s  may  be  written  as  (A  +  Z  )  sin  BAB  da. 
Hence,  the  ratio  of  the  cross-sectional  areas  may  be  expressed  as: 

(R  +  R  )2  sin  i p  dip  da 

y  =  ~ — h — — -  (143) 

(A  +  Z  )  sinflcostfi  dflda 
6  2  3 


[  12]  B.  van  der  Pol  and  H.  Bremmer,  "Further  Note  on  the  Propagation  of  Radio  Wave 
Over  a  Finitely  Conducting  Spherical  Earth,"  Phil.  Mag. ,  No.  182,  March  1939. 
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Or,  the  divergence  factor  Is  equal  to: 


D  =  {7 


VR3 

Ae+Z2 


f  sin  #  <ty 
sin  0cos  ^  d0 

u 


Expressing  d  d  in  terms  of  $  and  fi  permits  D  to  be  expressed  in  the  following  con¬ 
venient  manner.  Reference  may  be  made  to  the  original  works  for  the  details  involved. 

Ae(Rl  +  R2*  */  Sin  42  008  h 

D=  . . . = . —  . =  — . ======  (14! 

V  )  (Ae  +  Z2)R1  008  4 3  +  (Ae  +  Z1)R3  C0S  41  |  (Ae  +  Zl)(Ae  +  Z2)Sin  6 


Equation  (145)  is  a  function  of  the  variables  previously  expressed  as  desirable.  While 
it  may  be  slightly  unwieldy,  the  author  has  chosen  to  leave  it  in  its  present  form  so 
as  to  be  of  the  same  order  of  exactness  as  the  previously  derived  expressions  for  AR. 
In  general,  D  is  seen  to  have  its  smallest  value  when  the  reflected  energy  approaches 
a  grazing  angle  of  incidence.  As  $  decreases,  D  will  approach  unity. 

C.  The  Implementing  Expressions 

It  now  remains  for  the  functions  which  implement  the  previous  expression  to  be 
derived.  Referring  again  to  Figure  1 ,  the  tangent  line  designated  as  L  is  needed  to 
determine  when  the  functions  are  valid.  By  this  it  is  meant  that  when  Z  is  located 

mt 

within  the  interference  region  or  above  line-of-sight  of  the  transmitter,  then  the  pre¬ 
vious  equations  are  valid.  If  the  value  of  Z  as  combined  with  a  value  of  R  places  the 

«  Li 

receiver  or  target  out  of  sight  of  the  transmitter,  then  the  equations  are  not  valid.  Note 
that  a  line-of-sight  correspondence  between  the  two  is  also  a  degenerate  situation.  For 
a  given  value  of  transmitter  height,  the  distance  to  the  horizon  may  be  found  in  the 
following  manner.  From  Figure  1,  it  will  be  seen  that  a  radius  vector  drawn  from  the 
center  of  the  spherical  earth  intersects  the  tangent  ray  from  a  height  Z^  such  that  the 
angle  between  these  two  is  equal  to  ninety  degrees.  Hence,  the  vector  whose  magnitude 
equals  the  radius  of  the  earth,  the  tangent  ray,  and  the  vector  from  the  center  of  the 
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earth  to  the  transmitter  form  a  right  triangle.  Therefore,  the  distance  from  the  trans¬ 
mitter  to  the  horizon  can  be  expressed  in  the  following  manner: 


2  2 
+  A  )  -  A 
e  e 


(146) 


If  the  value  of  R^,  an  independent  variable,  is  greater  than  d^,  it  is  required  to  know 
at  what  height  ,  the  receiver  or  target,  must  be  in  order  that  it  falls  into  the  inter 
ference  region  or  where  the  equations  are  valid.  Consider  the  limiting  case  where  R 


is  at  some  point  tangent  to  the  surface  of  the  earth.  Call  the  distance  from  the  point 
of  tangency  to  the  projected  vector  from  the  center  of  the  earth  to  the  receiver  or 


target,  d: 


dl  =  R2-dh 


(147) 


A  radius  vector  from  the  earth '  s  center  to  the  tangent  point  is  again  at  right  angles  to 
d^  or  (R^  -  d  ) .  Hence,  the  intersection  of  the  vector  from  the  earth' s  center  to  the 
receiver  by  the  vector  d^  occurs  at  a  magnitude  of  the  target-directed  vector  equal  to 

d2: 

.  .2  ,2 

d2  ~  Ae  +  dl 


The  minimum  height  in  the  interference  region  is,  therefore: 


A  + 
e 


/ 


.2  £ 
A  +  d, 
e  1 


(148) 


In  order  to  specify  f(  <p  )  and  f(<?  )  as  a  function  of  R  ,  Z  ,  and  Z  ,  the  following 

1  4-  2  12 

analysis  may  be  carried  out  in  reference  to  Figure  7  which  is  taken  directly  from 
Figure  1.  The  analysis  Involves  first  the  derivation  of  6.  From  an  oblique  triangle 
consideration  define  a  variable  T  and  U  as  follows: 


T  = 


A  +  Z  +  A  +  R„  2A+R+Z 
e  2  e _ 3  e  3  2 


4 


(T  -  R3)(T  -  Ae)(T  -  Ag  -  Z2) 


(149) 


(150) 
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Hence,  Equations  (149)  and  (151)  serve  to  solve  for  6,  having  previously  determined 
R  .  The  following  relationships  follow  directly  from  Figure  1: 

u 


^3  =  *  -  6-0 

(152) 

iP4=ir  -  6-  flj 

(153) 

7T  7T 

—  a  =  —  -  a 

2  12  2 

(154) 

(155) 


The  value  of  y  in  Figure  1  is  determined  in  a  manner  similar  to  Equation  (151): 


where 


Hence  ^  =  y  -  \  4>2  \ 


y  =  2  tan 


-1 


(S  -  R2)(S  -  Rj) 


S(S 


*3> 


s  = 


+  R. 


2 


(156) 


0=2  tan 


-1 


(S  -  R2)(S  -  Rj) 


S(S  -  r3) 


(157) 


Within  the  preceding  framework  are  the  capabilities  of  analyzing  an  airborne-type 
transmitter  with  respect  to  its  pattern  of  coverage.  The  analysis  has  been  carried  out 
by  the  method  of  ray  analysis  of  the  electromagnetic  energy  leaving  the  transmitting 
antenna.  The  equations  are,  of  course,  subject  to  their  noted  limitations. 
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CHAPTER  V 

LIMITATIONS  OF  THE  ANALYSIS 

At  the  start  of  this  paper,  the  author  mentioned  that  all  equations  were  based 
upon  the  ability  to  represent  electromagnetic  waves  by  straight  lines  or  rays.  The 
author  also  indicated  that  this  assumption  may  not  always  be  exact  or  true.  In  the  pre¬ 
vious  analysis,  it  was  assumed  that  at  all  times  the  rays  appear  as  normals  to  the 
surfaces  of  constant  phase  of  the  electromagnetic  wavefronts.  However,  in  the  vicinity 
of  sharp  corners  and  in  regions  where  the  rays  exhibit  odd  properties  because  of  re¬ 
fraction.  such  as  at  a  cusp,  then  the  energy  can  no  longer  be  said  to  follow  the  rays  and 
geometrical  optics  fail  to  give  results  which  are  truly  meaningful.  At  this  point,  it  is 
required  to  analyze  the  energy  by  introducing  the  field  of  physical  optics.  Physical 
optics  is  not  the  theme  of  this  paper  but  the  following  analysis  should  indicate  to  the 
reader  where  the  previous  analysis  does  break  down  and  if  it  is  important  in  the  general 
area  of  pattern  coverage.  Also  covered  in  this  section  will  be  surface  roughness  and 
the  effect  of  the  index  of  refraction  of  the  atmosphere. 

A.  The  Equation  of  the  Ray 

Geometrical  optics  may  be  invoked  to  analyze  the  situation  where  the  atmosphere 
is  non -homogeneous.  In  a  homogeneous  atmosphere,  as  previously  assumed,  the  rays 
indeed  represent  the  wavefronts  and  their  direction  of  travel.  In  any  atmosphere, 
Maxwell 1  s  equations  apply,  and  hence  they  will  be  stated  for  a  desired  end  result: 


V.  D  =  p 

(158) 

V  .  B  =  0 

(159) 

(160) 

_  -  _  9D 

VxH  =  <TE  +  -aT 

(161) 

1=  mh 

(162) 
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D  =  €  E  (163) 

The  parameters  are  defined  as  follows:  D  is  the  electric  flux  density,  p  is  the  charge 
density  at  the  point  in  question,  B  is  the  magnetic  flux  density,  E  is  the  electric  field 
intensity,  H  is  the  magnetic  field  intensity,  p  is  the  permeability,  and  e  the  dielectric 
constant.  In  free  space  a  equals  zero  and  p  and  e  are  not  a  function  of  time,  and 
hence: 

V  x  E  =  -  p  yj-  (164) 

_  -  af 

VxH=£gr  (165) 

Then,  the  wave  equations  of  free  space  may  be  determined: 


Vx(VxE)  =  -  pVx 


3H 

at 


But: 

Therefore: 


In  a  similar  manner: 


V(V 


E)  -  V2e  : 


(VxH) 


V(V-  E)  =  0 


V2E  =  pc 


V2H  =  PC 


(166) 


(167) 


Equations  (166)  and  (167)  are  the  wave  equations  in  free  space.  If  we  assume  that  all  of 

i  u>t 

the  propagated  energy  can  be  represented  in  a  sinusoidal  form,  then  e  can  be  used 

to  represent  the  time  variations  of  all  fields  and  currents  involved.  The  second  time 

2 

derivative  can  hence  be  represented  by  -w  where  w  is  the  radian  frequency  of  the 
propagated  wave: 


_2-  2— 
V  e  =  peui  E 


(168) 
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or: 

V2H=M£W2H  (169) 

It  now  behooves  us  to  find  a  solution  to  these  equations.  In  general,  the  solution  may 
be  sought  In  the  following  form.  Either  Equation  (168)  or  (169)  may  be  considered. 
Considering  Equation  (169),  the  solution  Is: 

H  =  A  £  j  WS  (170) 

A  and  S  are  real  functions  of  position.  If  Equation  (170)  is  substituted  into  Equation 
(169)  and  the  real  and  imaginary  parts  are  equated  to  zero,  the  following  results: 

(VS)2-  —2  -  rj2  =  0  (171) 

Aw 

V2S  +  =  0  (172) 

where: 

2 

17  =  pe 

2 

If  to  is  quite  large,  then  Equation  (171)  may  be  written  as: 

(VS)2  =  rj2  (173) 

The  solution  to  Equation  (173)  is  simply  to  find  a  function  whose  gradient  of  the  magni- 

2 

tude,  but  not  the  direction,  is  given  as  ??  .  Hence,  according  to  the  type  of  wavefront 
to  be  represented,  a  surface  may  be  chosen  where  at  every  point  we  assign  a  constant 
value  of  S,  say  SQ.  The  gradient  of  S  is  perpendicular  to  S  at  each  point  in  question 
and  Vs  may  be  completely  defined.  To  find  another  surface  S  =  SQ  +  S,  perpendiculars 
to  the  defined  surface  may  be  constructed  of  lengths  equal  to  Vs  /  77  where  tj  is  the  value 
of  the  index  of  refraction  at  the  point  in  question.  The  traces  of  the  constructed  per¬ 
pendicular  vectors  then  define  the  desired  surface.  Hence,  this  is  the  method  of 
producing  surfaces  of  constant  S.  As  Vs  —  0  in  the  limit,  Equation  (173)  is  satisfied. 
The  surfaces  of  constant  S  represent  the  electromagnetic  wavefront .  The  perpendicu¬ 
lars  represent  the  rays  used  in  the  previous  derivations.  The  differential  equation  of 
the  ray  will  now  be  determined. 


50 


Consider  the  surface  of  Figure  S  to  represent  a  wavefront . 


In  this  figure,  t  is  the  normal  unit  vector  to  S  and  hence  represents  the  ray.  At  a 
point  T  1  on  the  ray  t  at  a  distance  S  +  A  S,  t  +  dt  represents  the  direction  of  the  ray 
tor  tne  new  surface  at  Xs  .  At  point  T'  ,  the  vector  representing  T'  is  t  ds/j];  as 
|  VS  |  =  t)  and  ds/rj  is  in  the  direction  of  Vs.  Hence  at  T  '  ,  t  +  dt  can  be  drawn  per¬ 
pendicular  to  S  =  S  +  A  s.  Let  u  and  v  be  unit  vectors  perpendicular  to  t.  By  definition 
v  =  t  x  u.  T' '  is  a  point  along  u  a  distance  dp  from  R.  T"  1  is  in  the  direction  of  t 
at  a  distance  ds/(t]  +  Vrj  •  udp)  where  7 )+  Vrj  •  udp  represents  the  new  value  of 
Tj  •  T'  and  T' ' '  can  both  be  proven  to  lie  in  the  same  wavefront.  As  tj  is  assumed  not 
to  change  in  the  v  direction,  it  can  also  be  proven  that  t  +  dt  is  normal  to  the  plane 
specified  by  T1  T' 1 '  and  v.  Hence: 


RT' "  =  dpu  + 


dS 

V 


Vn  •  udp 


t 


(174) 
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where: 


rjM  *J»I  I  I  _  p*Jl|  f  t 


-  ET' 


(175) 


as  u  x  v  =  t  and  t  x  v  =  -u 


hence: 


t  +  dt  = 


T'T"  '  xv 
|t't'"  I 


t  + 


Vn  '  u 
V 


dou 


(176) 


dt_ 

do- 


Vn  '  u 
V 


(u) 


(177) 


Equation  (177)  is  the  differential  equation  of  the  ray.  While  its  derivation  may  leave 

[  13] 

something  to  be  desired,  the  author  refers  the  reader  to  Kerr  ,  which  suggests  the 
foregoing  analysis.  Equation  (177)  gives  a  method  for  obtaining  a  ray  type  analysis  from 
a  direct  consideration  of  Maxwell '  s  equations.  Hence,  knowing  a  point  R  on  the  energy 
wavefront ,  a  point  T'  can  be  determined  by  choosing  a  value  of  do\ 

RT  =  td<r  (178) 

The  direction  of  the  ray  at  T  1  is  equal  to: 

t(T')  =  t(R)  +  —  (R)  da  (179) 

dcr 

As  dc  —0,  the  segments  so  obtained  approach  the  continuous  ray. 


B.  The  Effect  of  3 

It  is  now  interesting  to  look  at  the  ray  pattern  for  first  a  homogeneous  free-space 
atmosphere  and  then  the  stratified  atmosphere.  For  the  homogeneous  atmosphere, 
clearly  Vjj  is  equal  to  zero.  Hence  from  Equation  (177)  dt/dcr  is  equal  to  zero. 

Hence  from  Equation  (179)  the  rays  then  can  be  seen  to  be  indeed  straight  lines. 


[131  Kerr,  op.  cit.  ,  pp.  41-50. 
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However,  for  the  time,  assume  that  the  physical  earth  conditions  may  be  repre¬ 
sented  by  an  index  of  refraction  which  changes  only  in  the  height  direction.  If  the 
earth  and  its  atmosphere  are  represented  by  the  conventional  x,  y,  z  coordinate  system, 
then  tj  is  everywhere  only  in  the  z  direction.  It  then  follows  from  the  ray  equations 
that  for  this  somewhat  idealized  case  of  a  stratified  atmosphere  that  the  rays  repre¬ 
senting  the  correct  energy  wavefronts  are  curves  lying  in  planes  passing  through  the 
z  axis.  This  can  be  better  seen  if  the  absolute  value  of  Equation  (177)  is  obtained: 


dt 

dcr 


(u) 


If  |V rj  |  is  considered  vertical,  then  |Vtj  |  ■  u  =  | V rj  |  .  Also,  if  tj  is  considered 
approximately  unity  then: 


where  R  is  the  radius  of  curvature.  Hence^or  a  constant  gradient  of  rj  the  rays  are 
approximately  arcs  of  circles. 


C .  The  Modified  Index  and  Effective  Ae 

Under  these  assumptions,  the  straight-line  ray  analysis  is  seen  to  fall  apart. 
However,  it  is  possible  to  introduce  a  new  parameter  called  the  modified  index  of 
refraction  N  .  The  author  will  mainly  state  the  following  results  obtained  from  Kerr. 

A  modified  index  of  refraction  is  obtained  by  considering  Snell '  s  Law  and  the  previously 
derived  ray  equations: 

N  =  (  1  +  ^  )  TJ  *  rj  +  —  (180) 

as  rj  at  1 . 

Where  as  rj  generally  decreases  with  height  Z,  where  A  is  the  radius  of  the  earth, 
the  modified  index  N  will  be  seen  to  increase  with  height.  When  tj  decreases  linearly 
with  height: 


(181) 
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where  the  rays  are  represented  by  straighHine  segments  instead  of  segments  of  arcs. 
The  previous  is  obtained  by  a  consideration  of  Maxwell' s  equations,  the  ray  tracing 
equations,  and  Snell*  s  Law. 

If  dr)/dZ  (Z)  is  small  compared  to  t)  ,  the  value  for  free  space,  then: 


where: 


N  =  TJ 


0 


(182) 


J_  _  ±  An+  ± 

A  T).  dz  A 
e  0 


(183) 


A  represents  an  effective  radius  of  the  earth,  and  t}q  a  1 .  Measurements  have  shown 
that,  with  the  exception  of  a  few  hundred  feet  close  to  the  earth,  the  gradient  of  rj  is 
such  that  the  effective  radius  of  the  earth  is: 


=  4/3  (True  Radius)  (184) 

Hence,  through  the  utilization  of  an  effective  earth' s  radius,  a  linear  profile  for  the 
index  of  refraction  may  be  taken  into  account  in  the  analysis  of  the  interference  region. 
For  a  profile  other  than  that  of  the  linear  case,  the  preceding  becomes  a  poor  assump¬ 
tion.  However,  the  conditions  may  be  so  varied  from  one  physical  situation  to  another 
that  it  is  impractical  to  attempt  to  obtain  any  type  of  a  general  qualitative  analysis.  If 
the  profile  for  the  index  of  refraction  cannot  be  assumed  linear,  then  the  analysis 
would  fall  under  a  special  classification  such  as  duct  propagation. 


D.  Additional  Limitations  on  a  Ray  Analysis 

In  the  solution  of  the  radar  wave  equation,  the  ray  analysis  was  derived  under  the 
assumption  that: 


v2a 


Aw 


(185) 
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Through  the  analysis  of  this  inequality  and  Equation  (142)  it  can  be  shown  that  Equation 
(185)  Implies  the  restriction  that  the  index  of  refraction  must  not  change  appreciably 
in  a  distance  equal  to  that  of  the  wave  length  of  propagation.  This  is  the  reason  that 
as  the  frequency  increases  the  assumption  of  Equation  (185)  becomes  increasingly  more 
accurate.  For  a  ray  analysis  it  has  been  found  that  the  preceding  is  in  general  valid 
from  physical  measurements. 


An  analysis  of  Equation  (142)  also  imposes  the  restriction  on  the  ray  analysis 
that  the  fractional  change  in  the  spacing,between  the  bundle  of  rays  representing  the 
traveling  energy,  must  be  small  when  compared  to  unity  over  a  distance  of  one  wave 
length.  In  essence  it  must  be  remembered  that  a  single  ray  has  no  meaning}but  rather 
the  single  ray  is  used  to  indicate  the  composite  bundles  of  rays  which  in  turn  represent 
the  energy  wavefronts.  Hence,  in  the  analysis,  this  condition  is  violated  whenever  it 
assumes  that  a  focus  exists.  In  considering  rays  which  are  curved,  this  condition  is 
also  violated  when  the  rays  undergo  a  change  in  sign  of  their  curvature.  To  analyze 
the  conditions  imposed  by  Equation  (142)  it  becomes  necessary  to  resort  to  physical 
rather  than  geometric  optics. 


The  previous  ray  analysis  was  presented  with  the  intent  to  stress  that  problem 
areas  exist  where  a  ray  analysis  in  the  interference  region  may  lead  to  erroneous  re¬ 
sults  and  conclusions.  The  reader  should  hence  keep  these  areas  of  validity,  and  the 
assumptions  that  are  made  in  the  analysis,in  mind  when  using  the  theoretical  geometri¬ 
cal  equations  derived  for  the  interference  region.  The  author  did  not  analyze  in  detail 
these  situations  where  the  ray  analysis  leads  one  into  a  meaningless  procedure.  The 
many  different  situations  which  may  arise  should  be  given  their  own  separate  analysis 
utilizing  physical  optics  or  the  wave  equations  and  physical  optics.  For  the  remaining 
section  of  this  paper,  it  will  be  assumed  that  the  index  of  refraction  changes  linearly 
with  height  and  that  the  wavefronts  may  be  analyzed  by  the  method  of  straight  ray 
analysis, assuming  the  modified  radius  of  the  earth  to  take  into  account  the  physical 
curvature  of  the  rays  assuming  a  linear  profile  for  tj  . 
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E.  The  Hough  Reflector 

At  this  point,  it  might  be  beneficial  to  analyze  the  effect  of  the  reflective  Burface 
roughness.  It  has  been  previously  assumed  that  the  surfacecould  be  represented  by 
an  even  or  smooth  boundary  with  no  irregularities.  Consider  a  roughness  whose  height 
is  equal  to  h.  Let  ip  represent  the  angle  of  incidence.  Reference  may  be  made  to 
Figure  9. 


Figure  9.  Phase  Difference  between  Rays  Reflected  from  Two 
Levels 

The  phase  difference  between  the  two  rays  is  equal  to: 

A  R  =  2  h  sin  ip  (186) 

which  corresponds  to  a  phase  difference  of: 

KAR=  —  sin  ip  •  (187) 

Hence,  in  terms  of  this  phase  difference,  the  surface  roughness  may  be  specified. 

When  the  phase  difference  equals  it  ,  the  surface  is  effectively  its  roughest.  There  is 
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no  established  criterion  on  just  what  constitutes  a  rough  surface  but  the  following  may 
be  chosen: 

hsin^  <  (188) 

O 

If  h  sin  ip  <  ~  this  surface  may  be  assumed  smooth. 

O 

Needless  to  say,  this  is  only  a  crude  quantitative  analysis  of  surface  roughness. 

If  the  irregularity  is  of  the  order  of  a  wave  length,  then  the  problem  of  refraction 
should  be  considered.  Equation  (188)  does,  however,  emphasize  that  the  reflected  energy 
is  less  affected  by  a  given  height  of  roughness  as  ip  decreases  or  as  X  increases.  For 
a  surface  for  which  the  roughness  has  to  be  represented  as  a  statistical  distribution  of 
irregularities,  the  method  of  analysis  is  undeveloped  at  this  time.  Truly  a  ray  analysis 
fails  for  this  situation  and  at  present  electromagnetic  wave  theory  has  trouble  handling 
the  boundary  conditions  involved  for  a  statistical  distribution  of  surfaces.  Hence, an 
analysis  of  this  type  at  the  present  state-of-the-art  is  made  by  use  of  rough  and  simpli¬ 
fied  refraction  theory  or  from  direct  physical  observations  and  attempted  correlation. 
However,  very  little  may  be  said  at  the  present  time  in  terms  of  a  definite  correlation 
between  roughness  and  reflection. 

There  has  been  some  work  done  by  assuming  that  the  surface  roughness  can  be 
represented  by  a  series  of  ellipses  within  which  a  given  phase  retardation  for  reflection 
can  be  assumed.  The  basis  of  this  assumption  is  that  a  single  reflected  ray  is  in  essence 
representing  radiation  from  the  actually  illuminated  surface.  In  reality  a  large  surface 
is  illuminated  by  the  incident  energy.  Currents  at  the  surface  are  then  induced  which 
in  turn  cause  radiation  in  all  directions.  At  a  specific  point  in  space,  the  sum  of  the 
induced  radiation  components  yields  the  reflected  field.  The  area  and  dimension  of 
each  ellipse  can  be  derived  in  terms  of  the  wave  length  of  radiation  involved,such  that 
the  phase  of  the  radiation  from  each  ellipse  is  essentially  constant.  If  the  surface 
roughness  is  consequently  larger  than  the  calculated  ellipse,  then  an  analysis  can  be 
carried  out  by  assuming  essentially  multiple  reflection  points  in  a  ray  analysis  and 
making  calculations  based  upon  this  assumption.  The  usefulness  of  this  technique, in 
general,may  be  limited  and  hence  the  mathematical  derivations  have  been  excluded. 
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However,  the  reader  should  become  aware  of  the  surface  roughness  problem  In  a 
theoretical  analysis.  Even  if  the  surface  roughness  over  a  given  area  of  reflection  in 
a  physical  situation  could  be  specified,  which  is  many  times  difficult,  the  analysis 
techniques  are  not  well  defined.  A  qualitative  analysis  for  each  physical  situation 
involved  may  be  made  in  respect  to  the  preceding,  but  to  attempt  to  obtain  a  quantitative 
type  of  analysis  is  almost  hopeless  in  view  of  the  parameters  involved  and  the  state-of- 
the-art.  It  appears  that  a  theoretical  analysis  should  be  carried  out  by  assuming  a 
smooth  earth  and  then  modifying  the  results  qualitatively  should  the  roughness  phenome¬ 
non  appear  to  be  prominent  in  the  overall  analysis. 


F.  Antenna  Phase  Pattern  Effect 

It  might  be  interesting  at  this  point  to  consider  the  situation  where  the  energy 

radiated  from  the  transmitting  antenna  is  not  of  the  same  phase  in  the  direction  of  4> 

and  <p  .  From  the  previous  analysis  of  F  the  following  conclusion  can  be  immediately 

drawn.  If  f($  )  or  f (<p  )  are  complex,  then  they  can  be  represented  as  a  complex  number 
1  ^ 

in  the  expression  of  F  .  Hence,  the  analysis  as  carried  out  in  the  derivations  need  to  be 
only  slightly  modified.  The  effect  of  such  circumstances  would  be  difficult  to  analyze 
in  general.  However,  one  would  expect  that  the  shape  of  the  lobes  of  coverage  of  the 
transmitting  equipment  would  be  changed.  In  other  words,  the  area  of  maximum  coverage 
could  be  governed  in  this  fashion. 

The  feasibility  of  being  able  to  specify  a  phase  difference  from  the  different  parts 
or  sections  of  the  antenna  is  a  different  problem.  So  far  as  the  author  knows,  no 
commercial  producers  of  high-powered  antennas  specify  a  phase  difference  in  radiation 
from  various  sectors  of  their  antennas. 
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CHAPTER  VI 

A  DIGITAL  COMPUTER  ANALYSIS 


It  Is  appropriate  at  this  time  to  present  a  method  which  is  applicable  to  finding 
the  coverage  diagram  of  a  transmitter  located  over  an  assumed  smooth  surface  and  in 
an  atmosphere  which  has  a  linear  profile  for  the  index  of  refraction.  The  type  of  trans¬ 
mission  will  be  assumed  as  two  way  or  the  coverage  diagrams  obtained  will  be  that  for 
a  radar.  The  method  will  involve  the  solution  of  the  equations  of  Chapter  VI  (A)  by  a 
high-speed  digital  computer  —  specifically  the  IBM  704  or  709  or  7090.  However,  the 
resulting  computer  program  should  give  considerable  help  to  the  reader  in  applying  the 
techniques  to  any  digital  computer. 


A.  The  Necessary  Equations 

The  necessary  equations  to  completely  specify  the  problem  and  plot  the  coverage 
diagram  are  as  follows.  Reference  should  be  made  to  Figure  1. 

1.  Equation  (30)  gives  the  free-space  radar  range  equation: 


R 


0 


(30) 


where  G  is  the  maximum  gain  of  the  transmitter,  \  the  wave  length  of  the  transmitted 
waves,  Wt  the  peak  power  output  of  the  transmitter  in  watts,  the  minimum  value 
of  useful  power  to  the  receiver,  and  cr  the  target  area  in  square  meters. 

2.  The  following  equation  is  useful  to  give  the  maximum  possible  value  of  R^ 
under  ideal  conditions: 


R2  =  2R0 
zmax  u 


(189) 
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3.  An  empirical  solution  for  the  radius  of  the  earth  as  obtained  by  Massachusetts 
Institute  of  Technology  '  s  Lincoln  Laboratory  will  be  used.  This  equation  assumes  an 
ellipsoidal  earth: 

Ag  =  1852(3438.26815  +  577804  cos(2L)  +  0. 01222  cos(4L)  +  Hfg)  (190) 

where  L  is  the  latitude  in  degrees  of  the  area  covered  in  reference  to  the  geometric 
north  pole,  and  Hfs  is  the  height  above  sea  level  in  meters  of  the  area  involved.  A 
linear  profile  for  the  index  of  refraction  will  be  assumed  in  all  calculations,  hence  the 
effective  earth  radius  will  be  4/3  the  calculated  value.  The  region  of  reflection  will 
be  assumed  spherical  after  the  above  computation. 

4.  The  two  equations  which  indicate  the  boundary  of  the  interference  region 
will  be  employed.  These  equations  previously  appeared  as  Equations  (146)  and  (148). 

The  distance  from  the  point  to  the  horizon  is: 

V  y  <VA/-*e  <*«> 

If  Rg  is  greater  than  d^: 


J2min 


-A  + 
e 


/ 


*2  ,2 
A  +  d„ 
e  1 


dl  =  R2'dh 

5.  The  value  of  the  angle  formed  with  the  center  of  the  earth  as  a  function  of 
Z  ,  Z  ,  and  R  is  a  desired  parameter.  From  Equation  (125)  the  angle  results: 

1  a  Z 


9 1  =  2  tan 


-1  <8-VZ1*8-Ae-Z2> 


S(S  -  R2) 


2A  +  R„  +  Z  +  Z„ 
e  2  1  2 


(125) 
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6. 

and  R  . 

<5 


7. 


where: 


8. 


9. 

as: 


The  value  of  the  tangent  of  a  and  a  is  needed  to  find  the  distances  R 

1  i  1 

Equations  (117)  and  (118)  rewritten  appear  as: 


tana 


sin  ( 0 1  -  8) 


1  cos (^  -  8)  -R/(R+  zp 


(117) 


tana  = 


sin  9 


2  -  cos  8  +  R/(R  +  Z  ) 

it 

The  values  of  R  and  R  are  hence  given  by  Equations  (121)  and  (122): 

1  o 


(118) 


( [(z-i)Myz)/Ae]  ( 


) 


cosa. 


I  A 
)  6 


(121) 


Z  =  cosfS^  -  8) 


^  (cosfi  -  1)  +  (Z  cos  0)/A  | 

R3  |  cos  »2  ^  Ae 


(122) 


The  value  of  6  is  given  by  Equation  (151)  as: 


6=2  tan 


-1  ,  <T  -  R3)(T  ' 


T(T  -  Ag  -  Z2) 


T  = 


R3  +  2Ae+Z2 


(151) 


From  Equation  (151),  the  values  of  ip  ,  ip  ,  ip^  and  <t>  may  be  determined 

u  4  Z  £» 


=  *  -  6-9 


(152) 


=  ,r  _  §  -  $ 

4  1 


(153) 
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JT  V 

tli  =  —  -  a  -  —  -  a 
v2  2  12  2 


(154) 


**2  “  *  2  "  V 


(155) 


10.  The  value  of  </>^  is  needed  for  completely  specifying  the  transmitter  func¬ 
tions  and  appears  as  Equation  (157): 


2  tan 


_j  /  (S  -  R2)(S  -  Rj) 
S(S  -  r3) 


k 


s  = 


Rl  +  R2+R3 


(157) 


11.  The  complex  values  of  the  coefficient  of  reflection  are  given  by: 


r  = 


-  .  /I  2 

e  sin  ip  -  e  -  cos 


?  sin^  -cos2^ 


(134) 


rh 


sin  ^  -  /TTcos2  ^ 


Sin^2  +i/V 


2  , 
cos  ip 


where  (v)  and  (h)  designate  vertical  and  horizontal  polarization  respectively  and: 
r,  equals  -  j60A.ff>with  units  as  previously  defined. 

12.  The  divergence  of  energy  from  a  curved  surface  is  taken  into  account 

by  D: 


(135) 


D  = 


Ae(Rl  +  Hg)  s/sin^  coso^ 


\/{(Ae  +  Z2)R1C0S  ^3  +  (Ae  +  Zi>R3C0S  ^4  }  (Ae  +  Zl)(Ae  +  Z2)Sin  6 \ 


(145) 
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13.  From  Equations  (32)  and  (41),  the  propagation  factor  F  Is  given  by: 


F=  l[f(^)  +  f(02)  FDc'lY] 


(191) 


where: 


,  2  7T 

y  = 


(192) 


ar  =  ri  +  r3-r2?  0 


(193) 


-iy 

c  =  cosy  -  j  siny 


(194) 


14.  It  will  be  found  convenient  to  specify  a  quantity  K  as: 

£t 


K„  = 


.5 


2  R, 


2 

2R„ 


(195) 


From  Equation  (28),  the  power  received  by  the  radar  receiver  is  equal  to: 


W 


rl 


GV  o  F4 
(4*)3R* 


W, 


(196) 


with  units  as  previously  defined.  Utilizing  Equation  (30): 


K. 


(197) 
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Therefore: 


,  \V  „  2.2  4 

iK  f!4=  —  —  -G--~  ■  — 
1K2F|  W  4tt  2  4 

r  (4*)  R 


W  4  2.2 
t  F  G  X  o' 


W  ,^3d4 
r  (4u)  R2 


W 


rl 


W 


(198) 


Hence,  from  Equation  (198),  the  target  as  a  function  of  Z  ,  Z2>  and  R2  is  within  the 
pattern  of  coverage  of  the  transmitter  if: 


|KF|4  £  1.0  (199) 

A 

The  target  is  outside  the  area  of  coverage  of  the  receiver  if: 

,  ,4 

|K2F|  <  1.0  (200) 

Hence,  the  pattern  of  coverage  is  determined  by  assuming  a  fictitious  target  over  a 
suitable  grid  representing  the  interference  region  and  then  determining  at  each  point 
if  Equation  (199)  or  (200)  prevails. 

The  preceding  equations  have  been  chosen  because  they  are  advantageous  for  the 
purpose  of  plotting  a  coverage  diagram.  First,  any  desired  degree  of  accuracy  may  be 
maintained.  Secondly,  they  are  functions  of  the  parameters  used  in  a  coverage  analysis. 
Thirdly,  they  represent,  in  as  far  as  is  feasible,  the  actual  situation  involved.  True, 
some  situations  have  been  idealized.  However,  the  composite  relations  are  worthy  of 
a  detailed  analysis  because  they  give  an  insight  into  the  various  relations  involved.  In 
many  cases,  they  approach  reasonably  well  the  practical  situation.  This  is  especially 
true  for  an  area  over  the  ocean.  This  is  because  the  surface  of  the  sea  does  not  differ 
from  the  theoretical  as  greatly  as  a  terrain.  Certainly,  they  provide  at  their  worst, 
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basic  building  blocks  upon  which  the  coverage  analysis  may  rest.  This  is  becoming 
increasingly  Important  in  establishing  high  frequency  criteria  for  design. 

B.  General  Computation  Logic 

The  general  logic  used  to  solve  the  previous  equations  is  to  digitally  trace  the 
pattern  of  coverage  as  formed  in  space.  The  logic  first  assumes  some  starting 
value  of  range  R  and  target  height  Z  .  At  this  set  value  of  R  ,  the  value  of  Z  is 

Z  6  0  6 

reduced  until  the  boundary  condition  between  Equations  (199)  and  (200)  is  satisfied.  The 
programmer  must  Indicate  whether  this  first  point  is  on  the  top  or  bottom  of  an  assumed 
lobe-type  coverage  pattern.  The  stated  necessary  equations  must  be  solved  to  test  for 
each  point  in  space  considered.  When  for  a  set  value  of  R  the  correct  value  of  Z  is 

it  it 

determined,  then  a  tabulation  of  the  variables  is  made  and  R^  is  increased  or  decreased 
according  to  the  point  being  on  the  bottom  or  the  top  of  a  lobe.  The  maximums  and 
minimums  of  the  coverage  pattern  are  tested  by  A  R  changes  and  by  changes  in  the  sign 
of  the  difference  between  two  successive  power  calculations  by  Equation  (198).  Practi¬ 
cally,  the  logic  is  as  shown  in  Figure  X 0- 
cu 


Figure  10.  PIctorol  Logic  of  Pattern  Tracing 
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Much  time  has  been  devoted  to  applying  and  perfecting  the  rough  outline  of  the 
logic  previously  presented.  Its  many  obscure  points  require  a  thesis  for  elaboration. 
A  rough  outline  of  the  program  written  by  the  author  will  be  presented.  Reference  to 
the  particular  sections  of  the  Fortran  listing  should  be  made  to  pages  73  through  87. 
The  Fortran  logic  is  a  working  logic  and  has  been  used  on  an  IBM  7090.  Some  of  the 
results  may  be  found  in  the  appendices.  The  following  is  a  list  of  the  major  program 
symbols  in  the  terms  of  symbols  used  in  this  paper: 


AE 

= 

radius  of 

ARTHO 

= 

horizon  r 

Z2US 

= 

Z2 

RAN 

= 

R2 

ANG(2 ) 

= 

e 

ANG 

ei 

TAG1 

= 

TAN  «i 

TAG2 

= 

TAN  a 

RANG1 

= 

Ri 

RANG3 

- 

R3 

RANG13 

= 

ar 

APH02 

= 

“2 

APHOl 

= 

a 

1 

GMR 

= 

Re  <rj 

GMI 

= 

lm(r) 

DE 

= 

D 

SI2 

= 

*2 

SI3 

= 

*3 

SI4 

= 

*4 

TATA1 

= 

M 

TATA2 

= 

<p 

2 

ATA 

= 

6 

FELD1 

= 

F 

dh 
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CALL  (  )  Indicates  a  double  precision  routine 

DPREAD  -  Double-precision  read  of  Data  Cards 


DPADD  -  "  "  + 

DPMPY  "  "  x 

DPROOT  "  "  N f 

DPSUB  "  "  - 

DPDIV  "  "  t 

ASDP  -  "  "  SIN' 

DPSIN  -  "  "  SIN 

DPCOS  -  "  "  COS 


Referring  to  page  73,  Section  A  of  the  Fortran  listing  is  a  prologue  to  the  main  pro¬ 
gram.  It  assigns  storage  space  for  double -precision  quantities  which  are  used  in  the 
program.  Double-precision  operations  maintain  at  least  fourteen  significant  decimal 
digits.  These  are  required  at  points  previously  noted.  The  general  accuracy  of  a 
digital  computer  the  size  of  the  IBM  7090  is  seven  significant  decimal  digits.  Note 
that  all  sections  of  the  Fortran  listing  precede  the  symbol. 

Section  B  of  the  program  reads  in  data  stored  on  punched  cards.  Data  Card  1 
must  contain  in  order  of  their  appearance  in  fields: 

H  -  the  maximum  gain  of  the  transmitting  antenna. 

B  -  the  wave  length  of  propagated  wave  in  meters. 

D  -  the  peak  output  power  of  the  transmitter  in  watts. 

E  -  the  minimum  useful  signal  to  the  receiver  in  watts. 

F  -  the  effective  target  area  in  MKS  units. 

This  first  card  must  consist  of  five  data  fields;  each  field  is  ten  columns  wide.  The 
data  are  entered  in  the  floating  point  mode.  In  each  field,  the  following  must  hold: 

Column  1  -  sign  of  fraction. 

Column  2  -  optional  zero  punch. 

Column  3  -  the  decimal  point. 

Columns  4-9  -  the  quantity  expressed  as  a  decimal. 

Column  10  -  the  alpha-symbol  E. 
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Column  11  -  the  sign  of  the  decimal  exponent. 

Columns  12-13  -  the  power-of-ten  decimal  multiplier. 

Data  Cards  2  through  27  must  contain  the  value  of  the  antenna  pattern  stored  at  inter¬ 
vals  of  one  degree  in  elevation  from  +90  degrees  to  -90  degrees.  Starting  with  the 
first  card,  seven  values  of  the  antenna  pattern  function  are  stored  at  one  degree  incre¬ 
ments  starting  at  plus  ninety  degrees.  The  field  width  of  each  value  of  f (<p)  is  ten 
columns.  The  magnitude  of  f(  <p)  must  be  stored  as  a  three  digit  decimal  which  must 
occupy  columns  four  through  six  of  each  field.  The  remaining  columns  of  each  field 
are  utilized  as  explained  before.  Data  Card  28  must  contain  in  order  of  the  fields: 

Q  -  the  relative  dielectric  constant  of  the  reflecting  surface. 

R  -  the  conductivity  of  the  reflecting  surface. 

GAM  -  the  value  1. 0  if  horizontal  polarization  is  being  considered  and  the 
value  0.0  if  vertical  polarization  is  being  considered. 

AZ  -  the  azimuth  heading  for  which  the  180  degrees  values  of  i(<p)  are 
stored. 

TA  -  the  value  +1 . 0  to  test  lobe  logic  at  start  of  run  or  the  value  -1.0  if 
the  test  is  not  desired. 

ERD  -  the  value  of  production  block  printouts  desired  per  page. 

SLZ2US  -  the  slope  of  the  assumed  starting  lobe  as  referenced  to  rectangular 
coordinates. 

The  format  of  the  data  is  again  as  a  decimal  to  the  desired  power  of  ten  in  fields  ten 
columns  wide.  Data  Card  29  has  the  following  information: 

GMDO  -  has  the  value  of  -1.0  if  the  principal  root  is  to  be  used  in  determining 
T  .  The  value  +1.0  utilizes  the  program  as  written. 

STRAN  -  the  minimum  value  of  R^  allowed  in  program. 

RBCC  -  a  value  of  -  .  5  allows  for  the  grid  to  hunt  over  an  area  in  terms  of 

the  propagated  wave  length  equal  to  one  half  of  a  revolution  before  indi¬ 
cating  the  end  of  a  lobe.  A  value  greater  than  -  .  5  allows  less  area 
and  a  value  less  than  -  .  5  allows  for  a  greater  area. 

PRV  -  allows  for  monitor  printout  to  be  printed  with  tabulated  lobe  data. 
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Data  Card  30  consists  of  the  following  five  quantities: 

G  -  latitude  of  reflection  point  from  the  north  pole. 

ZAM  -  maximum  value  of  Z^  allowed  in  program. 

AIZ  -  the  minimum  value  of  allowed  in  program. 

RWHY  -  the  value  of  2.0  starts  the  program  over  a  lobe  tracing  logic  assumed 
as  the  bottom  of  a  lobe.  The  value  of  0. 0  assumes  the  starting  point 
as  the  top  of  a  lobe. 

EROE  -  the  value  of  4/3  modifies  Ag  as  assumed  by  a  linear  profile  of  r/ . 

The  next  six  data  cards  must  contain  double-precision  quantities  stored  one  per  card. 
They  appear  in  the  order  of  cards  as: 

AMD  -  the  maximum  allowable  difference  in  tan  and  tan 
relative  to  6. 

P  -  the  value  of  Z  . 

BET  -Az2- 

DEL  -AR  . 

w 

SRAN  -  starting  R  -AR  . 

SZ2US  -  starting  Z0  -  5  x  A  Z0 . 

The  format  of  these  cards  is  as  follows: 

Column  1  -  blank. 

Column  2  -  optional  zero  punch. 

Column  3  -  optional  decimal  point. 

Columns  4-20  -  sixteen  decimal  digits. 

Column  21  -  alpha-symbol  E. 

Column  22  -  the  sign  of  the  decimal  exponent. 

Columns  23-24  -  the  decimal  power  of  ten. 

The  statements  immediately  following  statement  15  calculate  R^,  2R^,  and  the  radius  of 
the  earth  by  Equations  (30),  (129),  and  (190)  respectively.  The  high  order  part  of  all 
double-precision  data  will  also  be  seen  to  be  set  equal  to  zero.  Statements  408  to  the 
end  of  Section  B  give  a  prologue  printout  of  the  program  parameters  on  tape  and/or 


paper. 
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Section  C  consists  of  pre-setting  the  lobe-tracing  logic. 

Starting  on  page  75,  Section  D  increments  Z  and  R  from  their  read-in  starting 

2  I* 


than  4,  Z  is  then  compared  to  Z  as  calculated  by  Equation  (148).  After  forcing 
“  2  2min 

Z  to  satisfy  this  condition,  program  stops  and  error  checks  are  tested.  Optional 
printout  routines  will  be  noted.  It  should  be  noted  that  in  Section  D,  statements  31  and 
32  increment  R2  in  the  proper  manner  when  this  section  is  branched  into  from  later 
sections  of  the  program.  Z2  is  likewise  adjusted  by  statements  39  and  40. 

A  double-precision  calculation  of  01  is  performed  in  Section  E.  Statements  478 
minus  1  through  statement  480  give  an  optional  monitor  of  the  critical  parameters  in¬ 
volved  in  this  calculation. 


Section  F,  covering  some  one  hundred  statements  through  statement  924,  serves 
to  find  two  values  of  6  which  may  then  be  used  in  the  straight  line  approximation  to 
make  a  equal  a  Section  F  starts  by  finding  a  value  of  $  by  the  flat  earth  approxi- 
mation.  This  appears  as  statements  830  through  831.  It  should  be  noted  that  throughout 
this  and  the  previous  section,  the  angles  themselves  have  not  been  calculated.  The  sine 
and  cosine  of  all  angles  are  used  directly.  This  procedure  was  used  to  increase  speed 
while  maintaining  the  desired  accuracy  throughout  the  program. 

After  having  found  a  value  of  6  by  the  flat  earth  approximation,  control  is  advanced 
some  eighty  statements  to  statement  130.  Statements  130  through  145  calculate  the  values 
of  the  tangents  of  and  c t  knowing  9  .  Note  that  this  is  essentially  a  manipulation  of 

Equations  (117)  and  (118).  Having  made  this  calculation,  control  is  transferred  back  to 
statement  832  and  through  to  statement  901 .  Monitor  action  will  be  noted  at  this  point. 

If  the  values  of  the  two  tangents  are  both  positive,  then  a  second  value  of  9  is  calculated 
by  the  incremented  flat  earth  approximation  or  Equation  (131)  slightly  modified.  This 
may  be  seen  in  statements  833  to  833  plus  3  on  page  78.  The  tangents  of  a  and  a 
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are  then  calculated  as  before  and  monitor  printout  will  be  noted  following  statement  845. 


At  any  time  during  the  preceding  two  approximations  should  the  angles  become 
such  that  the  tangents  of  a  or  a  become  negative,  then  the  following  logic  is  used. 

The  section  of  the  reflecting  earth  which  is  visible  to  both  the  transmitting  platform 
and  the  assumed  target  is  first  determined.  This  section  is  then  divided  into  four  equal 
sections.  Two  values  of  0  are  then  derived  which  causes  the  reflection  point  to  fall 
within  two  sectors  which  are  not  equal.  Having  two  valid  sets  of  data,  the  straight-line 
approximation  is  then  used  to  find  the  true  reflection  point.  If  the  previous  method  does 
not  yield  points  which  are  valid,  then  the  particular  point  in  space  under  consideration 
is  rejected  and  a  point  one  increment  in  height  is  considered.  This  forcing  procedure 
is  contained  in  the  statements  on  page  77. 

Section  G,  statements  148  to  421 ,  on  pages  78  and  79  calculates  the  new  value 
of  6  by  the  straight-line  approximation  such  that  tan  a  equals  tan  a  .  Optional  moni- 
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toring  printout  will  be  noted  as  statements  481  through  484  plus  1 .  The  validity  of  the 
calculations  are  again  attested  to  by  statements  495  and  936.  Statement  937  determines 
whether  the  two  tangents  are  equal  to  the  precision  desired.  If  they  are  not,  the  pro¬ 
cedure  is  repeated  start* ng  with  statement  148;  if  they  are,  the  program  continues. 
Statements  130  through  145  of  Section  F  are  used  in  the  calculation  of  the  tangents. 

Section  H  on  page  79  calculates  a  ,  a  ,  AR,  R  and  R  by  Equations  (193), 
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(121)  and  (122).  Optional  monitoring  is  again  provided. 

The  parameters  tp  ,  ip  ,  ip  ,  8  ,  a,  T,  and  <t>  are  calculated  in  Section  I.  The 
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calculations  are  straightforward  as  appearing  in  this  text.  The  polarization  of  T  is 
chosen  from  a  data  card. 

The  divergence  factor  D  is  calculated  by  Equation  (145)  in  Section  J. 

Section  K  on  page  81  converts  A  R  to  its  equivalent  phase  constant  by  Equation 
(192).  The  operations  of  Equation  (194)  are  also  performed. 

Section  L  reads  in  from  stored  data  the  values  of  f (<p  )  and  f(</>  )  and  interpolates 
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between  the  one  degree  stored  increments. 

Section  M  calculates  F  by  Equation  (191),  K  by  Equation  (195)  and  |k  f|  by 
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Equation  (198). 
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Within  Section  N  are  the  lobe-extremes  detecting  criteria.  This  section  covers 
branching  of  the  logic  at  lobe  maximums  and  minimums. 

Section  O,  starting  on  page  82  and  ending  on  page  84,  contains  the  lobe  tracing 
and  printout  logic  for  points  being  considered  as  nearer  the  top  of  a  lobe  than  to  the 
bottom.  This  section  decides  whether  the  value  of  Z  and  R  combined  comprise  a 
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valid  printout  point  and  if  not,  to  what  part  of  the  program  control  must  be  transferred. 
Z  and  R  are  increased  or  decreased  according  to  the  logical  decision.  This  section 
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makes  use  of  Equations  (199)  and  (200).  Statements  171  to  611  on  page  83  predict  the 
new  point  on  the  lobe  from  previous  information  of  valid  printout  points. 

Section  P  on  pages  84  to  86  is  a  duplicate  of  Section  O  for  points  on  the  bottom 
of  a  lobe. 

The  program  error  stop  printouts  and  formats  are  contained  in  Section  Q.  The 
section  also  contains  debugging  sequences  and  the  final  program  stop. 

The  fact  that  the  Fortran  listing  is  reduced  in  its  print  size  does  not  indicate  that 
it  is  to  be  taken  for  granted.  Within  its  framework  lies  the  answer  to  the  many  problem 
areas  that  arise  in  solving  the  equations  necessary  in  plotting  a  theoretical  coverage 
diagram.  Due  to  a  lack  of  space,  the  many  fine  points  and  problem  areas  have  not  been 
discussed  but  may  be  determined  from  the  Fortran  listing  by  the  individual  who  is  faced 
with  a  similar  problem.  It  should  be  noted  that  a  small  digital  computer  would  not  con¬ 
tain  the  storage  needed  by  the  logic  represented  in  the  previous  analysis.  The  program 
tested  maintains  an  overall  accuracy  of  at  least  three  decimal  digits  when  the  wave 
length  of  the  propagated  waves  is  equal  to  or  greater  than  one  tenth  of  a  meter. 
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Fortran  Listing 


DIMENSION  AMD(2),P(2),BET(2),DEL(2),AE(2),STl(2),SRAN(2) 

DIMENSION  ST2(2),ARTH0(2),ST3(2),Z2U3(2),RAN(2),FRED(2),C1(2) 
DIMENSION  Z2M(2),BB(2),C2(2),S(2),TM(2),TE(2),CB(2),C3(2),ANG(2) 
DIMENSION  SIl(2),Sl(2),CO(2),Z(2),ZB(2),DA(2),TAGl(2) 

DIMENSION  DB(2),TAG2(2),FRY(2),RAY(2),AVA(2),TRY(2) 

DIMENSION  CElM2(2),DLIM2(2),AFH0l(2),SL0C(2) 

DIMENSION  RANGl(2),APH02(2),RANG3(2),RANG13(2),A(l80),SZ2US(2) 
DIMENSION  VSANG(2),C0ANG(2),AEP(2),AEZ(2),SIANG(2) 

COMMON  SLOC 
SENSE  LIGHT  0 
PRINT  98 

98  FORMAT (1H1/61H  SENSE  SWITCH  0,1,0,1,0,1***PRESS  START  AND  CONTINU 
XE  IROGRAM) 

PAUSE 

IF  ACCUMULATOR  OVERFLOW  7,7 
7  IF  QUOTIENT  OVERFLOW  6,6 
6  IF  DIVIDE  CHECK  1,1 

1  READ  INFUT  TAPE  2,2,H,B,D,E,F 

READ  INPUT  TAPE  2,361, (A(j), J=l,l8o) 

READ  INPUT  TAPE  2,250,Q,R,GAM,AZ,TA,ERD,SLZ2US 
READ  INPUT  TAPE  2,628,GMDO,STRAN,RBCC,IRV 
READ  INPUT  TAPE  2,4,G,ZAM,AIZ,RWHY,ER0E 

READ  INFUT  TAPE  2,3,AMD(l),P(l),BET(l),DEL(l),SRAN(l),SZ2US(l) 
IF(ZAM)9,9,15 
9  READ  2,H,B,D,E,F 

2  FORMAT(5E13-6) 

READ  361,(A(J),J=1,1Q0) 

361  FORMAT (7E10. 3) 

READ  250,Q,R,GAM,AZ,TA,ERD,SLZ2US 
250  FORMAT (7E10. 3) 

READ  628,GMDO,STRAN,RBCC,PRV 
628  FORMAT (UE17- 9) 

READ  U,G,ZAM,AIZ,RWHY,ERQE 
4  FORMAT(E12 . 5 , 2E17 • 9 > 2E1.  5) 

READ  3,AMD(1),  P(l),BET(l),  DEL(l),  SRAN(l),  SZ2US(l) 

3  FORMAT(E24.9) 

15  0=4.0*3-14159265 
AB=(H*B)/C 

ac=(d/e)*(f/c) 

ad=sqrtf(ac) 

range=(sqrtf(ab) )*(sqrtf(ad) ) 

RANGM=2 . 0*RANGE 

AE(l) =1852.0* (3438. 2632-5. 778o4*(COSF(2.*G))-0.01222*(COSF(4.*G))) 
AE(2)=0.0 

CALL  DPADD(AE,AE,ST1) 

CALL  DPADD(ST1,P,ST2) 

CALL  DFMPY(ST2,P,ST1) 

CALL  DPR00T(ST1,ARTH0) 

sae=ae(i) 

ae(i)=ae(i)*eroe 

P(2)  =0.0 

AMD  ( 2 )  =0 . 0 
BET  ( 2 ) =0 . 0 
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DEL  (2)  *0.0 
NO-O 

SRAN  (2)*0.0 
SZ2US(2)*0.0 

IF  (SENSE  SWITCH  1)408,1(09 
1*08  FRINT  1*01 

1+01  FORMAT (55H  LOBE  PATTERN  PROPAGATION  TRACE) 

HUNT  4O2,AZ,ARTH0(l) 

1*02  FORMAT (17H0AZIMUTH  HEADING=E10.3,30H  LINE  OF  SIGHT  OF  TRANSMITTER* 
XE17.9) 

PRINT  1+07,  AE(l),  RANGE,  RAN0M,SZ2US(l) 

1*07  FORMAT (4HOAE=E17. 9, 7H  RANGE=E17 -9,7H  RANGM=E17.9,17H  STARTING  HEIG 
XHT=E17 • 9 ) 

PRINT  404,H,B,D,E,F 

404  FORMAT (10H0MAX  GAIN=E13.6,13H  WAVE  LENGTH=E1 3- 6,1011  PEAK  FWR=E13.6 
X,9H  MIN  SIG=E13.6,13H  TARGET  AREA=E13-6) 

FRINT  405,G,ZAM,Q,R 

405  FORMAT ( 24II0LATITUDE  FROM  NTH  POLE=E12.5,12H  MAX  TGT  HT=E17.9,19H  R 
XEL  DIEL  C0NSTANT=E10.3,l4H  CONDUCTIVm=E10.3) 

FRINT  4o6,P(l),BET(l),DEL(l),SRAN(l),AMD(l) 

1*06  FORMAT (12H0PLTF0RM  IIT=E17 .9,10H  DELT  ALT=E10.3,10H  DELT  RAN-E12.5, 
XI 4H  STARTING  RAN=E17.9,5H  AMD=E10.3/lHl) 

409  WRITE  OUTPUT  TAPE  7,401 

'.-/RITE  OUTPUT  TAPE  7, l*02,AZ,ARTH0(l) 

WRITE  OUTPUT  TAPE  7,407,AE(l),RANGE,RANGM,SZ2US(l) 

WRITE  OUTPUT  TAIE  7,404,H,B,D,E,F 
'./RITE  OUTPUT  TAPE  7, 405,G,ZAM, Q,R 

WRITE  OUTPUT  TAPE  7,4o6,P(l),BET(l),DEL(l),SRAN(l),AMD(l) 

'./RITE  OUTPUT  TAH2  15,401 

WRITE  OUTPUT  TAPE  15,402,AZ,ARTH0(l) 

WRITE  OUTPUT  TAPE  15,407,AE(l),RANGE,RANGM,SZ2US(l) 

'./RITE  OUTPUT  TAPE  15,404,H,B,D,E,F 
WRITE  OUTPUT  TAPE  15,405,G,ZAM,Q,R 

V/RITE  OUTPUT  TAPE  15,406,P(l),BET(l),DEL(l),SRAN(l),AMD(l) 

V/RITE  OUTPUT  TAPE  15,403 
5  CALL  DPADD(AE,P,AEP) 

CALL  DPDIV(AE,AEP,DA) 

C2(l)*2.0 

C2(2)=0.0 

RED=-30. 

KJ=-1 

C3(l)=1.0 

C3(2)=0.0 

GAY=0.0 

RAN(l)=SRAN(l) 

RAN(2)=SRAN(2) 

DRE=0.0 
R13P=0.0 
KC=0 
SBCC=0.0 
N0S=-1 
KLIEM*-1 
SSI=-1.0 
8  l/HY=RWHY 


YHW-1.0 

Cl(l)-4.0 

Cl(2)-0.0 

CALL  DPMFY(C1,BET,ST1) 

CALL  DPADD(9T1,SZ2U8,Z2U8) 

00  TO  30 

29  CALL  DPSUB(Z2US,BET,Z2US) 

30  IF(WH7-7HW)31,31,32 

31  CALL  DP8UB(RAN,DEL,RAN) 

00  TO  625 

32  CALL  DPADD(RAN,IXL,RAN) 

625  IF(RAH(1)-STRAN)156,33»33 

33  IP(RAN(l)-RAlBM)34,3Ml 

34  IP  (RAN-ARTHO)  35,10,36 

10  IP  (RAH(2)-AMH0(2))  35,35,36 

35  Z2M(l)«0.0 
zac(2)-o.o 
00  TO  37 

36  CALL  DP6UB(RAB,ARTH0.BB) 

CALL  DIMPY(AX,AX,STl) 

CALL  DPMFX(HB,HB,ST2) 

CALL  DPADD(8T1,3T2,8T3) 

CALL  DmOOT(9T3,8Tl) 

CALL  DPSUB(STl,AX,Za<) 

37  BCC-PHS 
FRED-  Z2US(l) 

FRED(2)-Z2US(2) 

SRC— 200 
AKA-0.0 

SEISE  LIGHT  1 
GO  TO  39 

38  SENSE  LIGHT  0 

39  CALL  DPADD(Z2US,BET,Z2U5) 

GO  TO  540 

40  CALL  DPSUB(Z2US,BET,Z2US) 

540  IP  (Z2US-ZEM)  460,541,42 

541  IF  (Z2US(2)-Z2M(2))  460,42,42 

460  CONTINUE 

462  IF (SENSE  SWITCH  3)463,464 

463  PRINT  6l4, Z2M(l),Z2US(l),RAN(l) 

614  FORMAT (5HOZ2M-E17. 9, 6H  Z2US-E17-9,5H  RAN-B17-9) 

464  WRIT*  OUTPUT  TAPE  7,6l4,Z2K(l),Z2US(l),RAN(l) 

00  TO  39 

41  IF  (SENSE  SWITCH  1)465,466 

465  PRINT  6l5,RAN(l) 

615  FORMAT (1H1/5H  RAN-E17-9,35H  RAH  IS  (HEATER  THAN  MAX  RANGE  STOP) 

466  WRITE  OUTPUT  TAPE  7,615, RAH(l) 

00  TO  156 

42  IF(Z2US(1)-ZAM)543, 543,43 

43  IF  (SENSE  SWITCH  1)467,468 

467  PRINT  6l6,Z2U8(l) 

616  FORMAT (1H1/6H  Z2U8-E17-9,37H  Z2US  IS  0REA1SR  THAN  MAX  TOT  HT  STOP) 

468  WRITE  OUTPUT  TAPE  7,6l6,Z2US(l) 
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00  TO  750 
543  CONTINUE 

CALL  DB(PY(ST1,ST1,ST1) 

CALL  DFMFY(RAN,RAN,ST2) 

CALL  DP6UB(8T2,ST1,ST1) 

CALL  DPADD(AE,Z2U8,AEZ) 

CALL  DTMFY(AEZ,C2,ST2) 

CALL  DFMPY(ST2,AEP,ST2) 

CALL  DPDIV(ST1,ST2,V3AB0) 

C  ,  C06(THETA(1))»1-VBAB3«C0AHj 

CALL  DPSUB(C3,VSAN0,CQAN0) 

C  SIN(THETA(1)^SOBT7(1.CQANG**2)-S«OT(1-COANO)(1+CQANO) 

CALL  DPADD(C3,CQAN0,ST1)  ' 

CALL  DB<PY(VSAN0,ST1,ST1) 

CALL  DHW0T(ST1,SIAW) 

ANO-ATANK (SIANO/CQANO ) 

C  DB«AE/(AE+Z2US ) 

CALL  DPDIV(AE,AEZ,nB) 

IP(SENSE  SWITCH  6)1+77,1*78 
1+78  IF(SENSE  SWITCH  5)1+79,1+80 
479  hunt  619,HAN(1),Z2US(1),AN0(1),AH02(1) 

619  FORMAT(5HORAN«E17.9,6H  Z2U8-E17.9,5H  ANO-E17.9,6h  AN02-B17.0) 

47?  £2  Sf 3 TAK  9) 

IND— 1 

IF(ssi)830,836,836 
836  ASSIGN  832  TO  K 
GO  TO  130 

830  GKA-(P(1)+Z2US(1))/ran(1) 

gkb«atahp(gka) 

GEC-1 . 57079632-GEB 
GED-Z2US ( 1 ) /OEC 

C  ,  .  SIN(THETA)-THETA 

si(i)-ged/az(i) 

Sl(2)-0.0 

831  ASSIGN  832  TO  K 
GO  TO  130 

83a  ASSIGN  147  TO  K 

IF(SENSE  LIGHT  3)835,845 
835  m(l)=TRY(l) 

FRY(2)-TRY(2) 

if (sense  switch  6)901,902 

902  IF(SENSE  SWITCH  5)903,904 

903  HUNT  632,TAGl(l),TAG2(l),m(l) 

632  FORMAT (52HO*******FLAT  EARTH  APFRCDtIMATIONtt»n"""««»iiii»Miiimi»iiii  6h 
XTAG1«E17.9,6h  TAG2-E17.9,5H  J7U.E17.9)  ’ 

H?HiT  631,SI(1),SI1(1),C0(1),Z(1),ZB(1) 

904  WRITE  OUTPUT  TAJS  7,632,TA0l(l),TA02(l),m(l) 

WRITE  OUTPUT  TAK  7,631, Sl(l),SIl(l),CO(l),Z(l),ZB(l) 

901  IF(TAG1)920,920,921  '  W'  V  ’’  K  1 

981  IF(TAD2)920, 920,923 


c 


r?  tangent  negative 


920  ir(lMD)927,925,925 

C  SIN(RAY)-SIN(ANG)C06(ANG-RAY)-SIN(ANC-RAY)C08(ANa) 

927  CALL  D»tFY(SIANO,AE,STl) 

CALL  DIMPY(CQA1»,ARTH0,ST2) 

CALL  DPSUB(STl,ST2,STl) 

C  IT  ARP  LINE  OP  SIGHT  OVERSHOOT  OF  TARGET*** . «hhhhhhhhhh>» . 

IF(ST1)940,  9ltO,94l 
940  SI-0.0 
SI(2)«0.0 
SI1-0.0 

sn(i)-o.o 

co-1.0 

C0(2)-0.0 
GO  TO  942 

9^1  CALL  DPDIV(STl,AEP,Sl) 

c  cos(ray)-cos(ang)cos(aic-ray)+sin(ang)sin(ang-ray) 

CALL  DFMPY(CQANG,AE,ST1) 

CALL  DBCPY(SIANO,ARTHO,ST2) 

CALL  DPADD(ST1,ST2,ST1) 

CAPP  DFDIV(STl,AEP,CO) 

C  AVA-SIN(Y)/3. 

C  »((SQRT(Z2US(2.*AE+Z20S))*C0-SI*AE)/(AE+Z2US)*3. 

C  IF  TGT  LINE  OF  SIGHT  OVERSHOOT  OF  ARP***************************** 
9^2  CALL  DPADD(AEZ,AE,STl) 

CALL  DIWPY(ST1,Z2US,ST1) 

CALL  DIHOOT  ( ST1 ,  ST1 ) 

CALL  DPDIV(ST1,AEZ,ST1) 

IF  (STl-SIAHG  )9**5 , 943 , 944 

943  IF(ST1(2)-SIAHG(2))945,944,944 

944  ST1-SIANG 
ST1(2).  SIANG(2) 

ST2.C0ANG 

ST2(2)-C0ANG(2) 

GO  TO  946 

945  CALL  DPDIV (AE,AEZ,ST2) 

946  CALL  DPMPY(ST1,C0,ST1) 

CALL  DIHPY(ST2,SI,ST2) 

CALL  DPSUB ( ST1 , ST2 ,  ST1 ) 

C  SINE  OF  THE  DIFFERENCE  DIVIDED  BY  . . . 

AVA-STl/3 . 

AVA(2)-0.0 
SENSE  LIGHT  3 
IF(AVA)925,925,926 
C  SI-SI+AVA 

928  CALL  DPADD(SI,AVA,SI) 

IND-IND+1 

GO  TO  831 

925  PRINT  930,RAN(1),Z2US(1) 

93O  FORMAT (36HCUNABIE  TOFIND  REFIECTION  POINT  RAN-E17-9,6H  Z2US-K17-9) 
933  WRITE  OUTPUT  TAPE  7,930,RAN(l),Z2US(l) 

GO  TO  38 
923  SI1(1)-SI(1) 
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SI1(2)-SI(2) 

i?( 180)929,926, 929 
929  IT(SSI)833,83M34 
834  SI(l)-(SIl(l)/lOOO.)+SIl(l) 

Sl(2)-(SIl(2)/lOOO. ) +8 11(2) 

GO  TO  831 

833  GEF«AE*((l.O/P)+(l.O/Z2U8) ) 

Sl(l)»(TRY(l)/(lOOO.O*GfcF) )+SIl(l) 

Sl(2)-0.0 
GO  TO  831 

845  IF(SKNEB  SWITCH  6)846,907 

907  IF  (SENSE  SWITCH  5)908,909 

908  PRINT  630,TAGl(l),TA02(l),FRY(l) 

630  FORMAT ( 52HO*HHHHHH>B(Cf(nigNTgJ  FLAT  KARTH  APPROXIMATION*** . *,6H 

XTAG1-K17-9,6H  TA02-E17-9.5H  FRY-K17-9)  . 

HUNT  631, 81(1), 811(1), C0(1),Z(1),ZB(1) 

631  FCRMAT(4H  SI-K17.9.5H  SH-U7-9>4h  CO-H17-9,3H  Z*S17-9,4e  28-E17-9  X) 

909  WRITS  OUTPUT  TAP*  7,630, TA0l(l),TAG2(l),m(l) 

WRITS  OUTPUT  TAPS  7,631, Sl(l),SIl(l),CO(l),Z(l),ZB(l) 

846  IF(TAGl)920,92O,924 
924  IF(TAG2)920,92O,l48 

C  C0-85PTF(1-8I**2) 

130  CALL  DWPY(SI,8I,ST1) 

CALL  DPSUB(C3,ST1,ST1) 

CALL  DPROOT(8T1,CO) 

C  ZB-8IANG*C0-SI*C0ANG 

CALL  DMPY(SIAHG,C0,ST1) 

CALL  DWPr(8I,C0A»,8T2) 

CALL  DFSUB(STl,ST2,m) 

C  Z-C0*CQANG+SI«3IANG 

CALL  DPMP3f(C0,CQANC,STl) 

CALL  DPMPr(3I,SIANB,3T2) 

CALL  DPADD(ST1,3T2,Z) 

C  TAG1-ZB/(Z- (AB/(AE+P) ) ) 

CALL  DPSUB(Z,0A,ST1) 

CALL  DPDIV(ZB,ST1,TAG1) 

C  TAG2»Sl/(C0-(AE/(AE+Z2US) ) ) 

CALL  DPSUB(C0,DB,ST1) 

CALL  DH)IV(SI,ST1,IAG2) 

NRC-NRC+1 

IF  (NHC)l45,l45,l63 
145  CALL  DPSUB(TAG1,TAG2,TRY) 

GO  TO  K, (147,832) 

147  CONTINUE 

AVA-IBY*( (81-811 )/(FRY-TRY) ) 

CALCULATE  NEW  SI 

148  CALL  DPSUB(SI,3I1,CE1M2) 

CALL  DPSUB(FRY,TRY,DL1M2) 

CALL  DPDIV(CE1M2,DL1M2,3T1) 

CALL  DFKPY(TRY,ST1,AVA) 

ST1-AMD*SI 

SI1-SI 

SI1(2)-SI(2) 


TOY-TRY 

TOY(2)-TRY(2) 

CALL  DPADD(AVA,SI.SI) 

48l  IF  (SENSE  SWITCH  6)495,1(62 

462  IF  (SENSE  SWITCH  5)463,464 

483  pribt  62o,si(i),sn(i),co(i),z(i),m(i) 

620  PORMAT(4h06I-S17*9,5H  SI1-K17.9,4h  CO-117-9, 3H  Z-E17-9,4h  TO-EI7.9  X) 

ranrr  62i,rayi  ,ava(i),toy(i),eaoi(i) 

621  PCFMAT(6h  RAY1-E17-9,5H  AVA-B17 •  9, 5H  TOY-117 . 9, 6H  TA01-E17.9) 

484  WRITE  OUTPUT  TAPE  7,620,Sl(l),SIl(l),C0(l),Z(l),ZB(l) 

WRITE  OUTPUT  TAPE  7,621,RAY1  ,AVA(l),TOY(l),TAOl(l) 

00  TO  495 

495  IP(TAQ1)920,920,936 

936  IF(TAG2)920,920,937 

937  IP  (ABSF  (AVA)  -  ST1)  151,151,130 

C  COS  (APH01 )  -1 . 0/SQRT  ( 1 . 0+TAQl**2 . )  ^t1«^^*^««hhhhh,#«hhhhhhhhhh, 

151  CALL  DFMPY(TAG1,TAG1,ST1) 

CALL  DPADD(ST1,C3,ST2) 

CALL  DFROOT ( ST2 , ST1 ) 

CALL  DPDIV(C3,ST1,ST1) 

C  C08(APH02)-1.0/SORT(1.0+TA02**2. ) «ST3 ' **»«»' » ** »*« * ** » *•  ****** ********* 
CALL  OPMPY(TA02,TA02,ST2) 

CALL  DPADD(ST2,C3,ST2) 

CALL  DPR00T(ST2,ST2) 

CALL  DPDiy(C3,ST2,ST3) 

C  RAN01-AK ( ( (Z-’l  .0)>(P*Z)/AE)/C06(AFH01)  )**  ****** . . . «»*♦« 

CALL  DPMPY(Z,P,ST2) 

CALL  DPDIY (ST2, AE,ST2) 

CALL  DFSUB(Z,C3,ST4) 

CALL  DPADD(ST4,ST2,ST2) 

CALL  DPDIY(ST2,ST1,ST2) 

CALL  DPHPY(ST2,AK,RANG1) 

C  RANG3-AK((  (C0S(RAY)-l.O)+(Z2*COS(RAY)  )/AE  )/C0S(APHO2)  )««**#****#****  | 

CALL  DFMPY(Z2US,C0,ST2) 

CALL  DPDIV(ST2,AK,ST2) 

CALL  DPSUB(C0,C3,3T4) 

CALL  DPADD(ST2,8T4,3T2) 

CALL  DFDIV(ST2,ST3,ST2) 

CALL  DPHPY(ST2,AE,RAI»3) 

C  . . .  »«««*«  ««  ••  *•  *««  «  ********** . 

CALL  DPADD(RAN01,RAN03,Sr4) 

CALL  DFSUB(ST4,RAN,RANG13) 

AFHO-ATAN?( (TA31+TA02)/2 . 0 ) 

APH01-ATAN7 (TAC1 ) 

APH02-ATANP(TA02) 

PHS-RAH013/B 

IF (SENSE  SWITCH  2)465,486 

486  IF (SENSE  SWITCH  5)467,488 

487  PRINT  622,TAGl(l),ATO0l(l),TA02(l),APn02(l) 

622  FORMAT(6HOTAG1-E17.9,7H  AFH01-E17.9,21H  TA02-E17-9, 

X7H  APH02-E17.9) 

438  WRITE  OUTPUT  TAPE  7,622, TA01(l),APHOl(l),TA02(l),APH02(l) 

485  ATA-3.14159265-APH0  1 


80 


SSI-+1.0 

512- 1.57079632-AJHO 
RAn-ATANF(SIl/CO) 

513- 3 • 14159265-AIA-RAYl 

Sl4»3 . IU159265-AIA- (AMO (l)-RAYl ) 

SA- (RANOl ( 1 ) +RAN03 (1 )+RAN (1 ) ) /2 . 0 
SB-(SA-RAN(1 ) )*(SA-RAN0l(l) ) 

SC-SA*(SA-RAM03(1)) 

SQAC-SQRTF(SB/SC) 

TATA2—  ( 1 . 57079632-814) 

TATAI-2 . 0# (ATANF (SQAC ) ) -ABSF (TATA2 ) 

BR6-60.0*B*R 

C2AEH0-C03F  (SI2)«COSF  (SI2 ) 

GA-Q-C2APH0 
OB— ABSF(BR6) 

gc-(ga«oa)+(gb«gb) 

GD-SQRTF(aC) 

BETA2-ATANF  (OB/OA) 

BETA1-BETA2/  2  ■  0 

gdr-sqrtf(gd) 

CHETAl-COSF(BBTAl) 

SBETAl-SIHF(BKTAl) 

GE»CBETA1*0DR 

GF«SBETA1*GDR 

DZ-0.5 

C  CURVED  AMD  FLAT  EARTH  DIFF«(RAN**2)/2.0*SAE  RAMOE-RAN*COS(TATAl) 
SUZ2US- (RAN**2 ) /(2 . 0*SAE ) 

PZ2US-Z2US-3UZ2US 

IF(0AM-DZ)152.152,153 

152  GG-Q*SIHF(SI2) 

GH— BR6*SHff(SI2) 

GI-GG-GE 

GJ-GH-GF 

GK-GG+GE 

GL-GH+GF 

gm»(gi*gi)+(gj*gj) 

gn-sqrtf(gm) 

BETA3-ATANF(GJ/0I) 

go«(gk*gk)+(gl*gl) 

GP-SqpTF(GO) 

betaWtanf(gl/ck) 

gq-gn/gp 

GR=BETA3-BETAl* 

GMVH=GMDO* (GQ*C  OGF(GR) ) 

GMVX^}Q*SINF(GR) 

GO  TO  155 

153  GS«SINF(SI2)-GE 
GT— GF 

0U-3INF(SI2)+GE 
GV«(GS*GS)  +  (GT«GT) 

GW-SQRTF(GV) 

BETA5-ATAMF(GT/GS) 

beta6»atanf (gf/gu) 


GX-  { GU*GU )  +  ( GF*GF ) 

GY-SQRTF(QX) 

OZ-GW/GY 
HA-BETA5-BETA6 
GMHR^JMDO*  (GZ*COSF  (HA) ) 
GMHI-GZ*SIKF(HA) 

155  HB«SINF(APHO)*COSF(AFHO) 

HC-SQRTF(HB) 

HD-RANOl (1 ) +RANG3 ( 1 ) 

HE«AE(l)*HD*HC 

HF- ( AE ( 1 ) +Z2US ( 1 ) ) *S  Uff  ( ANG ( 1 ) ) 
HC-(AE(l)+P(l))*HF 

HI- ( AE ( 1 ) +Z2US ( 1 ) ) *RAMG1 ( 1 ) *COSF ( S 13 ) 

HJ-(AE(1)+P(1))*RAN03(1)«C06F(SI4) 

HK-HJ+HI 

HL-HK*HG 

HM-SQPTF(HL) 

de-he/hm 

DELT-{RAHG13 (1)/B )*6 . 26318530 
DELH4-MODF  ( BELT , 6 . 26 318530) 
CSD-C06F(DELTM) 

SND-S INF ( DEL1M ) 

DZ-0.50 

IF(GAM-DZ)3U9, 3*+9, 350 
3^9  GMR-GMVR 
GMI-GMVI 
GO  TO  352 
350  GMR-GMER 
GMI-GMHI 
352  PB-GMR+CSD 
PC«GMR*SND 
PIWJMI*CSD 
H5-GMI*SHD 
PF-re+PE 
PG-PD-PC 

AITA2— ( (TATA2*57  •  2957T95  )-90 . 0) 
AITA1— ( (TATA1*57  •  2957795  )-90 . 0 ) 
J-XIHTF(AITA2) 

AL0T2«A(j) 

J-J+l 

FL0A2-FL0ATF ( J ) 

AL0T21-A(J) 

J-XINTF(AITAl) 

ALOTI-A(J) 

J-J+l 

FL0A1-F10ATF ( J ) 

ALOTll-A(j) 

ADDA2- ( AL0T21-AL0T2 )» (FL0A2-AITA2 )/l . 0 

ADDA1- (ALOT11-AL0T1 ) * ( FL0A1-AITA1 )/l . 0 

F0T1-AL0T11-ADDA1 

F0T2-AL0T21-ADDA2 

PR-F0T2*DE*PF 

PI-F0T2*I**PG 


FQEl-FOTl+fR 

FCE2-PI 

FOE-  (FCE1*FQE1 ) + (FCE2*FCB2  ) 

FELDl-SQPTT(FCE) 

AYK2.((O.5)*RAN0M)/RAN(l) 

SFWR-FWRl 
PAY1»AYK2*FELD1 
PWR1 -PAY1*PAY1*PAY1*PAY1 
IF(TA)115,U5,44 
44  if(why-yhw)68, 68,69 

69  IF(KC-2)68, 68, 70 

70  IF(FWR-1. 0)71, 115, 115 

71  IF(SFWR-FWB1)68,  68, 175 

66  if(bcc-(phs+rbcc))175,175,H5 
175  KJ=0 

IF(WHY-YHW)47,47,46 
1+6  IF(RANG13-R13P)l60,l60,48 

46  R13P*«RANG13 
CO  TO  75 

47  KC-0 
SBCC-SBCC+1.0 
GO  TO  50 

115  IF(WHY-YHW)51,51,?6 

50  WHY=2.0 

CALL  DPADD(RAN,DEL,RAN) 

Z2US(l)-FRE0+BET-(FZ2US*2. ) 

IF(SENSE  SWITCH  1)475,474 

474  IF(SENSE  SWITCH  5)29,476 

475  PRINT  6l8,RAN(l),Z2US(l),FRED(l),WHY,NRC 

6l8  FORMAT (31H0L0BE  CHANGE  EXIT  4  RAN»E17-9,6H  Z2US-E17.9,6H  FR 

XED=E17-9,5H  WHY=E10.3,5H  NRC-13) 

476  WRITE  OUTPUT  TAK  7,6l8,RAN(l),Z2US(l),FRED(l),WHY,HRC 
GO  TO  29 

51  IF(PWRl-1.0)6o,6l,6l 

60  IF(BET(l)-AKA)64,444,444 
444  KC-KC+1 

IF (SENSE  SWITCH  4)101,102 

102  IF  (SENSE  SWITCH  5)103,104 

103  PRINT  110,FWRl,RAN(l),Z2US(l),FBS,DE,GMR,GMI,FOTl,FOT2,NRC,IND 
110  FORMAT (17H0EXIT  10  10  FWR1«E9.2,5H  RAN—  3FF8 . 3 , 4H  Z2-OFF10.3,5H  FH 

XS-F8.2,4H  DE-F4.2,4H  GM-2F6.3,5H  F0T=2F5.2,5H  NRC-13,5H  IND-12) 

104  WRITE  OUTPUT  TAP®  7,110,PWRl,RAN(l),Z2US(l),PHS,DE,GMR,GMI,FOTl,FO 
XT2,NRC,IND 

GO  TO  40 

101  IF (SENSE  SWITCH  2)40,447 
447  WRITE  OUTPUT  TAPE  7,610 
6l0  FORMAT (llHCEXIT  10  10 ) 

IF (SENSE  SWITCH  3)446,449 
446  PT=1.0 
PRINT  610 
GO  TO  448 
64  TA-1.0 

IF(NOS)84,84,85 


84  SBCC-INTF(FHS) 

NOS.l 

85  NO-NO+1 
KC-0 

IF(SENSE  SWITCH  1)450,451 

450  PRINT  6o6,NO,PWRl,RAN(l),Z2US(l),PZ2US,PHS 
DRE-DHB+1.0 

IF (DRE-ERD) 451,860,860 
860  PRINT  607 
180  DRE-0.0 

451  WRITE  OUTPUT  TAPS  15,410, NO, PVRl,RAN(l),Z2US(l). PZ2US,PHS 
IF(SENSE  SWITCH  2)16,17 

16  IF(SENSE  SWITCH  4)18,19 

18  IF(FRV)19,22,22 

17  WRITE  OUTPUT  TAPE  7,6o6,N0,FWRl,RAN(l),Z2US(l),PZ2US 

WRITE  OUTPUT  TAPE  7,6O2,RAN0l(l),RA*)3(l),BA!C13(l),rH8,NRC,IHD 
WRITE  OUTPUT  TAPE  7,603,F0Tl,TAIAl,F0T2,TATA2,nOAl,FIOA2 
WRITE  OUTPUT  TAPE  7, 6o4,DK, OMR, OKI, 812,813,814 
GO  TO  22 

19  WRITE  OUTPUT  TATE  7,23,NO,TAOl(l),AFHOl(l),FOTl,FOT2,DK,GMR,GMI,NR 
XC,IND 

23  FORMAT (4H  NO-16, 6h  TAG1-E10.3,7H  AFH01-E10.3.6H  FOT1-E10.3,6H  F0T2 
X-E10.3,4H  DE-F4.2,5H  OMR-F6.4,5H  OMI-F6.4,5H  NRC-13,5H  IND-12) 

22  IF(KJ)171,172,173 

171  SFZ2US-Z2US 
KJ-1 

Z2US-Z2US+SLZ2US 

167  RED-REDfl.O 

IF (RED-ERD ) 31 , 6ll , 6ll 
173  FZ2US-Z2US-SFZ2US 
SFZ2US-Z2US 
Z2US-Z2US+FZ2US 
GO  TO  167 

172  SFZ2US-Z2US 
Z2US-Z2US-FZ2US 
KJ-1 

GO  TO  167 

6ll  WRITE  OUTPUT  TAIE  15,607 
WRITE  OUTPUT  TAPE  15,403 

403  FORMAT (12H  PRINT  NO. ,20H  RELATIVE  POWER, 17H  RANGE  IN  K 
XM,22H  HEIGHT  IN  METERS,  21H  PLOTTING  HEIGHT, 24H  MODU 
XLO  PHASE  SHIFl/lHO) 

RED— 30.0 
GO  TO  31 

6l  AKA-1 .0+BET(l) 

65  IF  (SENSE  LIGHT  1)66,67 

66  FRED-Z2US 

67  KC-KC+1 

IF(SENSE  SWITCH  4)27,28 
28  IF(SENSE  SWITCH  5)52,53 

52  PRINT  56,FWH1,RAN(1),Z2US(1),PHS,DE,GMR,0MI,F0T1,F0T2,NRC,IND 
58  FORMAT (17HCEXIT  11  U  PWR1-E9-2,5H  RAN— 3PF8.3,4H  Z2-0PF10.3,5H  PH 
XS-F8.2,4H  DE-F4.2,4H  GM-2F6.3,5H  FOT»2F5.2,5H  NRC-13,5H  IND-12) 


53  WRITE  OUTPUT  TAPS  7,58,WR1,RAN(1),Z2US(1),FHS,DE,OMR,OMI,FOI1,FOT 

x2,nrc,ind 

00  TO  3fl 

27  IF (SENSE  SWITCH  2)38.455 
455  WRITE  OUTPUT  TAP*  7,612 
612  FORMAT (11H0EX  IT  11  11) 

IF(SENSE  SWITCH  3)^54,457 
454  PT-1.0 
PRINT  612 
00  TO  456 

75  WHY-O. 0 

CALL  DFSUB(RAN,DEL,RAN) 

Z2USL-Z2US(l) 

Z2US-FRED+BET- (FZ2US*2 . ) 

IF  (SENSE  SWITCH  1)472,473 

472  PRINT  6l7,RAN(l),Z2USL, FRED (l), WHY, NRC,Z2U8(l) 

617  FORMAT (31HOLOBE  CHANGE  EXIT  14  RAN-B12.4,7H  Z2USL-E12.4,6H  F 

XRED-E12.4,5H  WHY-E9-2,5H  NRC-13,6H  Z2US-E12.4) 

473  WRITE  OUTPUT  TAPE  7,617, RAN(l),Z2USL,FRED(l), WHY, NRC,Z2US(l) 

GO  TO  30 

76  IF(PWR1-1. 0)78, 79,79 

78  IF(BET(l)-AKA)8l,92,92 

79  IF(KLIEM)87,87,86 
87  KC-0 

IF(Z2US(l)  +  (-BET-AIZ-Z2M) )82,80,80 

80  AKA-1 . O+BET ( 1 ) 

KC=KC+1 

IF(SENSE  SWITCH  4)105,106 

106  IF (SENSE  SWITCH  5)l07,108 

107  PRINT  111,FWR1,RAN(1),Z2US(1),PHS,DE,GMR,GMI,F0T1,F0T2,NRC,IND 
111  FORMAT (14HCEXIT  10  PWR1-E9.2,5H  RAN—  3PF8.3,4H  Z2-0PF10.3,5H  PHS-F 

X8.2.4H  DE-F4.2,4h  GK-2F6.3,5H  FOT-2F5-2,5H  NRC-13,5H  IND-12) 

108  WRITE  OUTPUT  TAPE  7,111, PVRl,RAR(l),Z2US(l),FHS,DE,GMR,GMI,F0Tl,F0 
XT2,NRC, IND 

GO  TO  40 

105  IF (SENSE  SWITCH  2)40,421 

421  IF (SENSE  SWITCH  3)422,423 

422  IRINT  601 
FT—1.0 

601  FORMAT (8HQEXIT  10) 

448  HUNT  605,PWR1,RAN(1),Z2US(1) 

605  F0RMAT(6H  FWR1-E17.9,6H  RAN«E17.9,6h  Z2US-E17-9) 

HUNT  602,RANGl(l),RANG3(l),RAl»13(l),HiS,NRC,IND 

602  FORMAT (7H  RANG1-E17-9,8H  RAN03«E17-9,9H  RAJKJ123«E17.9,5H  PHS-F8.2 
X,5H  NRC-13,5H  IND-12) 

HUNT  603,FOT1,TAIA1,FOT2,TATA2,FLOA1,FLOA2 

603  FORMAT (6H  F0T1-E12 . 5 , 7H  TATA1-E17.9,6H  F0T2-E12. 5,7H  TATA2-E17.9-6 
XH  J(1)=E10.3,6H  J(2)-E10.3) 

PRINT  604,DE,GMR,GMI,8I2,SI3,Sl4 

604  FORMAT (4H  DE-E17.9.5H  GMR-E12.5,5H  GMI-E12.5,5H  SI2«El6.8,5H  SI3-E 
Xl6.8,5H  SI4-E16.8) 

420  IF (PT) 423, 423, 449 

423  WRITE  OUTPUT  TAPE  7,601 


1*49  WRITE  OUTPUT  TAPE  7,605,WRl,RAH(l),Z2U8(l) 

WRITE  OUTPUT  TAPE  7,602,RA*Jl(l),RA»3(l),RA»13(l),PH8,HRC,IHD 
WRITS  OUTPUT  TAPE  7,6O3,FOTl,TAIAl,FOre,TAIA2,FL0lAl,FLQA2 
WRITE  OUTPUT  TAPE  7, 604,nE,GMR,0MI,SI2, 313,314 
GO  TO  40 

81  IF(SBCC-(PHS+0.5))83,S6,86 
86  Z2US-Z2US+BET 

KLHM-1 
GO  TO  540 
83  TA-1.0 
NO-NO+1 
KLHK--1 
KC-0 

82  IF(SKNSE  SWITCH  1)430,431 

430  HURT  6o6,N0,PWRl,RAH(l),Z2l»(l),PZ2U8,FHS 

606  FORMAT  (lflOl6,l6H  RELATIVE  POWKR-EI7.9.8H  RANGE—  3PS8 . 2, 11H  XM  HE  I 
XGHT-0PF8.2,9H  M  PLOT-78. 2,15H  M  PHASE  SHIFT-F13.7) 

ERE-ERE+1.0 

IF(nRE-KRD)431,459,459 

459  HURT  607 

607  FORMAT  (lHl) 

181  DRK-0.0 

431  WRITE  OUTPUT  TAPE  15,410, HO, PVRl,RAH(l),Z2UB(l),PZ2US,FHS 
4lO  FORMAT ( I10,E24 .9, -3FT13 . 2, 0FF20 . 2,F21 . 2,F25 -7) 

IF (SENSE  SWITCH  2)20,21 

20  IF (SENSE  SWITCH  4)24,25 

24  IF(PRV)25, 26,26 

21  WRITE  OUTPUT  TAPE  7,6o6,N0,FWRl,RAN(l),Z2US(l),PZ2US 

WRITE  OUTPUT  TAPE  7,602,RAHDl(l),RANG3(l),RANG13(l),PB3,NRC, HID 
WRITE  OUTPUT  TAPS  7,603,FOT1,TATA1,FOT2,TAIA2,FLOA1,FLQA2 
WRITE  OUTPUT  TAPE  7,6O4,0E,OMR,OMI,SI2,SI3,Sl4 
GO  TO  26 

25  WRITE  OUTPUT  TAPE  7,23,HO,TAOl(l),APHOl(l),70Tl,FOT2,DE,GMR,GMI,NR 

xc,ntD 

26  IF(KJ)i68,i69,170 

168  SFZ2US-Z2US 
KJ-1 

166  RED-RED+1.0 

IF (RED-ERD ) 32, 60S , 608 
170  FZ2US-Z2US-SFZ2US 
SFZ2US-Z2US 
Z2US-Z2UB+FZ2US 
00  TO  166 

169  SFZ2US-Z2US 
Z2US-Z2US-7Z2UB 
KJ-1 

00  TO  166 

608  WRITE  OUTPUT  TAPE  15,607 
WRITE  OUTPUT  TAPE  15,403 
RED— 30.0 

00  TO  32 
92  KLIEM—  1 

IF(SENSE  LIGHT  1)93,94 


93  FRED-Z2US 
91*  KC4CC+1 

IF(SKNSE  SWITCH  4)54,55 

55  IF(SENSE  SWITCH  5)56,57 

56  PRINT  59,FWR1,RAN(1),Z2U8(1),FHS,DE,0KR,GHI,F0T1,F0T2,KRC,IHD 

59  FORMAT (14HCBX  IT  11  FWR1-H9.2,5H  RAN— 3PF8 . 3, 4H  Z2-0P710.3, 5H  PHS-F 
X8.2,4H  nE-F4.2,4H  GM-2!F6.3,5H  FOT-2F5.2,5H  NRC-13,5H  IND-12) 

57  WRITE  OUTPUT  TAPE  7,59,WRl,RAN(l),Z2US(l),FHS,DE,OMR,GMI,FOTl,FOT 
X2,NRC, IND 

GO  TO  38 

54  IF(SENSE  SWITCH  2)38,458 
458  IF (SENSE  SWITCH  3)432,433 

432  PRINT  609 
PT— 1.0 

609  FORMAT (8H0EXIT  11 ) 

456  PRINT  605,FWR1,RAN(1),Z2US(1) 

PRINT  602,RANG1(1),RANG3(1),RANG13(1),PHS,NRC,DID 
HUNT  603 , F0T1 , TATA1 ,  F0T2,  TATA2,  FL0A1 ,  FLQA2 
PRINT  604,DE, OMR, GMI,SI2,SI3, Slit 

IF(PT)433,433,457 

433  WRITE  OUTPUT  TAPE  7,609 

457  WRITE  OUTPUT  TAPE  7,605,FWRl,RAN(l),Z2US(l) 

WRITE  OUTPUT  TAP*  7,602,RANGl(l),RANG3(l),RANG13(l),PHS,NRC,IND 
WRITE  OUTPUT  TAPE  7,603,FOT1,TATA1,IW2,TATA2,FLOA1,FLOA2 
WRITE  OUTPUT  TAPE  7,604,DE,0MR,GMI,SI2,SI3,Sl4 
GO  TO  38 

160  PRINT  l81t 

WRITE  OUTPUT  TAPS  7,l84 

184  FORMAT  (1H1/13H  LOWER  RAN013) 

ASSIGN  75  TO  K1 

161  PRINT  183,RAN(1),Z2US(1),RANG13(1),PU3P 

WRITE  OUTPUT  TAPE  7,l83,RAN(l),Z2US(l),RANG13(l),R13P 
183  FORMAT  (5H  RAN-  E17.9,6H  Z2US-  B17.9,8H  RANG13-  E17.9,6H  R13P-E1 
17-9  /lHO/53HO»»«»»»»»»»»»^»«»«»»«8IMUIAIKD  PROGRAM  STOP********** 
2/ltOHO  PROGRAMED  PDUMP  NCW  GOING  ON  COMJUNm/lHO/lHO) 
L-XL0CF(SL0C)-99 
CALL  PCUMP(SLOC(L),SLOC(1),3) 

PRINT  185 

WRITE  OUTPUT  TAPE  7,185 
END  FILE  7 
END  FILE  15 

185  FORMAT (91H0FR0GRAMED  PDUMP  CaMPIETED»»»»»««EXIT  . . . 

LEND  OF  PR0GRAM»«»lt»»»»»«FRED  LOWELL) 

CALL  EXIT 

162  PAUSE 

GO  TO  XI, (75,145) 

163  PRINT  102 


WRITE  OUTPUT  TAPE  7,182 
182  FORMAT  (1H1/26H  200  ANGIES  WITHOUT  OUTPUT) 
NRC-- 200 


ASSIGN  145  TO  K1 
156  GO  TO  l6l 

624  FORMAT (1H1/29H  END  OF  PROGRAM  F  LOWELL  RAN-E17-9,6h  Z2US-E17.9) 


750  mure  624,ran(i),z2us(i) 

END  FILE  7 
END  FILE  15 
CALL  EXIT 

END(l, 1,0, 0,1, 0,0, 0,0, 0,0, 0,0, 0,0) 


750  PRINT  624,RAN(1),Z2US(1) 
END  FILE  7 


89 


CHAPTER  VII 
CONCLUSION 


A  method  for  a  rigorous  theoretical  analysis  of  the  pattern  of  coverage  of  trans¬ 
mitting  equipment  in  the  earth1  s  atmosphere  requires  equations  which  are  cumbersome 
to  solve.  The  formulation  of  A  R  by  the  transcendental  approach  of  the  author  yields 
results  which  may  be  sought  to  a  high  degree  of  accuracy  with  a  minimum  of  work.  The 
overall  concepts  of  reflection,  divergence,  and  the  index  of  refraction  of  the  atmosphere 
may  be  taken  into  account.  The  analysis  has  been  shown  to  be  only  feasibly  solved  by 
a  digital  computer. 

Further  work  should  be  directed  to  one  of  two  areas.  The  first  is  that  of  obtaining 
an  index  of  refraction  other  than  that  which  yields  the  4/3  earth.  It  appears  that  no 
constant  value  can  be  assumed  throughout  the  region  of  entire  space.  Any  variances 
from  the  standard  4/3  earth  noticeably  causes  deviation  in  any  pattern  of  coverage.  The 
second  is  through  the  utilization  of  a  normalizing  procedure  of  the  equations.  This 
would  noticeably  reduce  the  work  required  in  performing  the  calculations. 
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APPENDICES 


APPENDIX  A 


The  following  two  graphs  should  supplement  Chapter  III  and 
Chapter V  of  this  report.  Figure  A-I  shows  the  dependence  of  AR  in 
^seconds  as  a  function  of  the  range  R  .  The  transmitter  and  target 
heights  are  maintained  at  constant  values.  Figure  A-l  was  obtained  by 
the  exact  equations. 

Figure  A-2  shows  the  dependence  of  A  R,  as  a  function  of  R  , 

a 

on  the  method  used  for  calculation.  The  exact  and  cubic  formulations  of 
AR  are  compared.  All  calculations  of  this  appendix  were  made  by  N.I. 
Durlach,  A.  M.  Carpenter,  and  M.A.  Herlin  of  Group  45  at  Lincoln 
Laboratory. 


A-l 


AR  (/a  SECONDS) 


Figure  A-l 

AR  in  /i seconds  as  a  function  of  Rj  as  calculated  by  the  exact  expressions  for 
a  standard  and  4/3  earth.  The  height  of  the  transmitter  is  40,000  feet  and  the 
height  of  the  receiver  or  target  is  60,000  feet. 
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Figure  A-2 

A  R  in  /^seconds  as  a  function  of  Rj  as  calculated  for  a  4/3  earth  by  the 
exact  equations  and  the  equations  represented  in  the  section  of  the  cubic 
r^.  The  height  of  the  transmitter  is  20,000  feet  and  the  height  of  the  tar¬ 
get  or  receiver  is  60,000  feet. 
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APPENDIX  B 


The  following  figure  and  tables  were  obtained  by  assuming  typical  radar  param¬ 
eters  as  follow.  The  maximum  gain  of  the  antenna  is  251,  the  wave  length  of  the 

g 

propagated  waves  is  0. 705  meters,  the  peak  power  output  of  the  transmitter  is  2  x  10 

*•15 

watts,  the  minimum  useful  signal  is  4  x  10  watts,  the  target  at  Z  has  an  effective 

A 

cross-sectional  area  of  10  square  meters,  the  latitude  from  the  north  pole  is  0.747 
radians,  the  relative  dielectric  constant  of  the  reflecting  surface  is  80,  the  conduc¬ 
tivity  is  4.  3  mhos  per  meter.  The  radius  of  the  earth  as  calculated  by  the  computer 
from  Equation  (190)  is  6, 366, 878.23  meters.  This  yields  an  effective  radius  of 
8,488,958-70  meters.  Horizontal  polarization  is  assumed. 

Utilizing  the  data  and  a  typical  antenna  pattern,  Figure  B-l  Is  sketched  from 
values  calculated  by  the  computer  logic  previously  listed.  The  computer  traced  out 
the  noted  lobes  of  Figure  B-l  with  A  R  equal  to  2 , 000  meters  and  A  Z  equal  to  2 . 5 

u 

meters.  Z^  equals  4,572  meters.  Figure  B-l  shows  the  coverage  of  the  second, 
third,  and  fourth  lobes  obtained  when  a  4  /  3  earth  was  assumed.  The  first  lobe  is 
omitted  because  it  is  very  close  to  the  surface  of  the  earth.  The  two  lobes  which  are 
larger  were  obtained  using  a  geometric  earth.  The  value  assumed  for  r\  is  hence  seen 
to  be  very  critical. 

Table  I  and  Table  II  list  some  of  the  major  factors  discussed  in  this  paper 
calculated  at  various  points  on  the  lobes  obtained  from  a  4/3  earth  consideration. 
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Figure  B-l 

Pattern  of  coverage  of  a  radar  with  typical  properties  stated  on  page  B-l .  Geo¬ 
metrical  and  4/3  earth  coverage  is  noted. 
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Table  I 

» 


Typical  intermediate  parameters  necessary  for  the  calculation  of  a  coverage 
diagram. 


L  x  103 

Z 

|k2f|4 

R  x  103 

~TT 

R3  x  10  j 

AR 

2 

2 

Eq.  (198) 

Eq.  (121) 

Eq.  (122) 

Eq.  (193) 

220 

55.5 

.347 

218 

9.79 

1.23 

240  ; 

85.5 

.584 

233 

14.2  i 

.907 

266  | 

135 

.504  j 

242 

23.0 

.855 

286  | 

225 

.944 

248 

39.4 

.962 

310 

390 

.980 

251 

58.2 

1.057 

286 

262 

.958 

244 

39.1 

1.22 

266 

182 

.442 

239 

26.8 

1.26 

246 

122 

.643 

229 

16.0  1 

1.21 

226 

ITS 

t— 

00 

.256 

1  216 

9.28  j 

| 

1.23 

206 

62.5 

.769 

198 

5.36 

1.34 

216 

84.0 

.303 

208 

9.25 

1.63 

236 

119 

.355 

222 

15.4 

1.59 

256 

174 

.344 

233 

24.6 

1.56 

276 

264 

,  .858 

239 

38.4 

1.64 

296 

389 

.938 

243 

54.0 

1.69 

296 

418 

.951 

241 

53.1 

1.88 

276 

298 

.965 

236 

38.6 

1.91 

256 

213 

.506 

228 

26.3 

2.02 

236 

153 

i 

.217 

217 

17.2 

2.13 

216 

j  113 

.324 

204 

10.7 

2.19 

206 

104 

.281 

198 

8.91 

2.23 

226 

141 

.285 

212 

14.5 

2.21 
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Table  II 


Typical  intermediate  parameters  necessary  for  the  calculation  of  a  coverage 
diagram. 


R2xl03 

Z2 

tan  a  ,  a 

1  2 

Eq.  (117) 

Eq.  (118) 

Re,rh> 

Eq.  (135) 

D 

Eq.  (145) 

I  m(r  )  X  10  3 

a 

Eq.  (135) 

226 

55.5 

123 

-.999 

.916 

+  .616 

246 

85.5 

170 

-.999 

.816 

+  .463 

266 

135 

218 

-.999 

.687 

+  .388 

286 

225 

264 

-.999 

.574 

+  .297 

310 

390 

299 

-.999 

.478 

+.257 

286 

262 

241 

-.999 

.579 

+  .318 

266 

182 

196 

-.999 

.685 

+.  391 

246 

122 

157 

-.999 

.798 

+  .499 

226 

87.5 

119 

-.999 

.890 

+  .654 

206 

62.5 

91.4 

-.998 

.947 

+.  855 

216 

84.0 

101 

-.998 

.916 

+.  759 

236 

119 

132 

-.999 

.837 

+.592 

256 

174 

167 

-.999 

.735 

+.470 

276 

264 

199 

-.999 

.633 

+  .394 

296 

389 

227 

-.999 

.550 

+.  344 

296 

418 

218 

-.999 

.554 

+.359 

276 

298 

187 

-.999 

.634 

+.  412 

256 

213 

156 

-.999 

.730 

+  .500 

236 

153 

126 

-.999 

.825 

+.621 

216 

113 

98.4 

-.998 

.903 

+  .784 

206 

104 

86.6 

-.998 

.923 

+  .  900 

226 

141 

109 

-.998 

.861 

* .  711 
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